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Abstract 

We study the electromagnetic form factor of the pion and kaons at low-energies 
with the use of Chiral Perturbation Theory. The analysis is performed within the 
three flavour framework and at next-to-next-to-leading order. We explain carefully 
all the relevant consistency checks on the expressions, present full analytical results 
for the pion form factor and describe all the assumptions in the analysis. From the 
phenomenological point of view we make use of our expression and the available 
data to obtain the charge radius of the pion obtaining (r 2 )^ = (0.452 ± 0.013)fm 2 , 
as well as the low-energy constant Lg(m p ) = (5.93 ± 0.43) x 10 -3 . We also obtain 
experimental values for 3 combinations of 0(p e ) constants. 
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1 Introduction 



The theory of the strong interaction is a gauge theory of quarks and gluons, Quantum 
Chromo Dynamics (QCD). At high energies and/or short- distances these degrees of freedom 
and a perturbation expansion in the gauge coupling constant work well and provide a lot 
of experimental tests of it. These provide the basis for our belief that QCD is indeed the 
theory of the strong interaction. 

At low energies and long distances, it is more difficult to extract experimental con- 
sequences out of QCD. An approach that has had a lot of successes is the use of chiral 
symmetry and its constraints as originally done using current algebra. This relies on 
the fact that QCD has a SU(nj)^x SU(nj)^ chiral symmetry in the limit of n^ massless 
flavours. This symmetry is spontaneously broken by nonperturbative QCD dynamics to 
the diagonal vector subgroup SU(3)y. 

Chiral Perturbation Theory (ChPT) is the effective field theory method to use this 
property at low energies. It takes into account the singularities associated with the Gold- 
stone Boson degrees of freedom caused by the spontaneous breakdown of chiral symmetry 
and parametrizes all the remaining freedom allowed by the chiral Ward identities in low 
energy constants (LECs). The LECs are the freedom in the parts of the amplitudes that 
depend analytically on the masses and momenta. The expansion is then ordered in terms 
of momenta, quark masses and external fields. Recent lectures introducing this area are 
given in Ref. JIJ. ChPT in its modern incarnation was founded by Weinberg || after 
earlier work by Pagels and collaborators [[J. Gasser and Leutwyler extended and system- 
atized Weinberg's work in a series of papers which caused the revival of chiral methods 
0, |H . The use of the external field method and dimensional regularization allowed a fully 
chiral invariant treatment throughout, obviously independent of any parametrization of 
the Goldstone Boson fields. 

The expansion is ordered in terms of momenta, meson masses and external fields. We 
use here the standard ChPT counting where the quark mass, scalar and pseudoscalar 
external fields are counted as two powers of momenta. Vector and axial-vector external 
currents count as one power of momentum. The Lagrangian can be ordered as 

£ 2 + £ 4 + C 6 + ■ ■ ■ 

10 2 90 94 

£ 2 +J2 Li o\+ ^2 H t o\ °i + J2 C * • -• ■ (LI) 

1=1 i=l 1=1 i=91 

The index i in £j stands for the chiral power. The precise form of £2 and £4 is given in 
Sect. [| while Cq can be found in ||. 

This expansion can be done both for the case of two or three light flavours, the number 
of terms given here correspond to the three light flavour case. The lowest order, 0(p 2 ), in 
this expansion corresponds to tree level diagrams with vertices from £ 2 , the next-to-leading 
order, NLO or C»(p 4 ), to one-loop diagrams with vertices from £ 2 or tree level diagrams 
with one vertex from £ 4 and the rest from £ 2 . The next-to-next-to-leading order, NNLO 
or 0(p 6 ), has two-loop diagrams, one-loop diagrams with one vertex from £4 and tree 



^effective 
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level diagrams with one vertex from L% or two vertices from £4 and all other vertices from 
£ 2 - The loop diagrams take all singularities due to the Goldstone Bosons correctly into 
account. The singularities are the real predictions of ChPT while the other effects from 
QCD are in the values of the LECs. When comparing ChPT and QCD with experiment, it 
is thus important to distinguish between the two types of contributions. In QCD one has 
to use a scheme where the regularization is consistent with chiral symmetry. The remnant 
of the precise definition at short distances in QCD reflects itself in the precise value of B , 
defined below, and the high-energy constants Hi, H 2 and Cj, % = 91, . . . , 94. 

The present situation is such that at lowest and next-to-leading order, ChPT is quite 
predictive since we can with relatively few assumptions determine all parameters and es- 
tablish relations between various observables. In the two-flavour case this program has 
been to a large extent carried to NNLO as well. In the three light flavour case this is 
still in progress. The present manuscript is one more step in the calculation of mesonic 
processes to NNLO. We present here the full NNLO calculation in three flavour ChPT of 
the pion and kaon electromagnetic form factors. In particular this allows us to determine 
one more low-energy parameter U 9 to this order. Additional assumptions on the 0(p 6 ) 
constants are at present still necessary and we will discuss some of these below as well as 
experimental determinations of three of them. 

The determination of the L\ is important since they are the real consequences of QCD at 
low energies, distinguishing it from other theories with the same chiral symmetry pattern. 
They are also needed to improve the predictability of ChPT. The previous fits of L\ were 
done at NLO where the dependence is at most linear and the correlations can thus be 
undone by using appropriate linear combinations of observables. In practice we still have 
large correlations present in L\, L r 2 and L\. At NNLO the dependence on the L[ can become 
quadratic and it is easier to perform the fits using MINUIT, taking correlations into account 
directly. This what we have done in our earlier work [|7|, [8|, ||, [H]] for L\, i = 1, . . . ,8. We 
now add L r 9 to this list. What we found earlier is that the added correlations at NNLO are 
present but not very strong, but that the dependence on the Zweig suppressed constants 
L4 and Lg can become substantial. We have nothing new to add on this question at present 
but will return to it in future work. Work at NLO exists, see [IT] and references therein. 



Earlier work on the electromagnetic form factors in ChPT is the NLO work of [12, 13], the 



calculation of the nonanalytic parts by [|14j] and the two-flavour calculation at NNLO [IS]. 



Some results for vector form factors at NNLO in three flavours also exist [If:, 17, 18 



Here we first give a few definitions of ChPT to fix our notation in Section |2|. In Sect. 
!T] we define the vector form factors. The known NLO results from [12, 13] are rederived 
and presented in Sect. |3.2| . The main new analytical part of this work, the vector form 
factors in the three flavour case in the isospin limit, is presented in Sect. |3.3| and compared 



as much as possible with the known results at this order [jlq, O, 18]. The data and our 



ChPT based fits are discussed in Section |j. The results for L\ to NNLO in ChPT accuracy 
as well as experimental determinations of 3 combinations of 0(p 6 ) constants are in Section 
|5[. This includes resonance estimates of the 0(p 6 ) constants needed in Sect. ^TT]as well as a 
comparison with predictions for L r 9 in Sect. |5.4| . In Sect. ^|we summarize our conclusions. 
We present some of the lengthier expressions for the pion form factor and the integrals 



3 



we have used in addition to those we discussed in J?| in the appendices. 



2 Some definitions 

The expressions for the first two terms in the expansion of the Lagrangian are given ex- 
plicitly by (Fq refers to the pseudoscalar decay constant) 

£2 = ^{(D^D^U) + ( X ] U + X U ] )} , (2.1) 

and 

£4 = L^D^D^U) 2 + L 2 (D^D U U){D^D V U) 

+L z {D^D i JJD v U^D v U) + U{D^D^J){^U + X U ] ) 
+L 5 (D^D II U( X j U + rf X )) + L 6 ( X j U + X U ] ) 2 
+L 7 ( X ] U - X U ] ) 2 + LstfUrfU + XU ] XU ] ) 
-iL 9 (F^ u D^UD u U j + F^D»U ] D U U) 

+L l0 (U^F^UF L n + HtiF^F*" + F^ u F Lflu ) + H 2 ( X j X ) , (2.2) 

while the next-to-next-to- leading order is a rather cumbersome expression |J . The special 
unitary matrix U contains the Goldstone boson fields 



U = exp I ) , M 



( >° + > vr+ K+ \ 



(2.3) 



The formalism we use is the external field method of H] with s, p, 1^ and r M matrix val- 
ued scalar, pseudo-scalar, left-handed and right handed vector external fields respectively. 
These show up in 

X = 2B (s + tp) , (2.4) 

in the covariant derivative 

DyP = d„U - ir^U + iUlp , (2.5) 

and in the field strength tensor 

= W» ~ d Ar), ~ * [i(r) M , Z(r)„] . (2.6) 
For our purpose it suffices to set 

\ / 2/3 \ 

(2.7) 

with e the absolute value of the electron charge and the classical photon field. 

Even if in Eq. ( |2.2| ) there is no explicit reference to the high energy regime, the theory 
depends on it via the values of the LECs. 
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3 The form factors: analytical results 
3.1 Definition 

The most general structure for the on-shell pseudoscalar-pseudoscalar-vector Green func- 
tion is dictated by Lorentz invariance. With the additional use of charge conjugation and 
electromagnetic gauge invariance one can parametrize the pion and kaon electromagnetic 
matrix elements as 

<7r + (g)b>+(p)> = (q, + P,)F v (t), 
(K+(q)\j„\K+(p)) = (q»+p»)F« + (t), 

(K°( q )\j,\K°(p)) = (q,+p,)Ff(t), (3.1) 

with t = {q — p) 2 - The current refers to the electromagnetic current of the light flavours 

2 1 - 

3n = 3 («7m m ) - 3 { d ln d + «7^s) • (3-2) 

The quantities Fy, Fy° and Fy + will be referred to hereafter as the vector form factors or 
simply the form factors. They are also defined in the crossed channel {0\jn\M a (p)M b (— q)). 

Notice that we explicitly neglect here the part of the electromagnetic current due to 
the heavier quarks, 

= I + h„t) - I (M) • (3-3) 

Its effects are suppressed by Zweig's rule and by the inverse of the heavy mass. 
These form factors are near t = often described by the charge radius, 



{r )v ~ b dt v 



(3.4) 



t=Q 



(3.5) 



and the quadratic term in the expansion, 

V " 2 dt^ v 
Electromagnetic gauge-invariance imposes the constraints 

i^(0) = if + (0) = l, F*°(0) = 0. (3.6) 

There is no corresponding form factor for 7r° and the rj. They vanish due to charge conju- 
gation. 
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Figure 1: One- loop diagrams contributing to the form factors. The square stands for the 
insertion of the vector external current. The cross indicates an 0(p A ) vertex from £4. 



3.2 Leading and next-to-leading order 

The leading order expressions for the vector form factors follow from the tree level diagrams 
with the Lagrangian ( |2.1| ) . The diagram is depicted in Fig. |I](a) . The result is rather simple 
and the form factors are constant, just those for a point-like particle, and always have the 
value (B1I), 



(t) 



0. 



(3.7) 



The simplicity of Eq. ( [3.7D is of great help when computing next-to-next-to-leading cor- 
rections. It will be sufficient to use next-to-leading expressions only when rewriting quark 
masses and chiral limit decay constants in terms of the physical meson masses and decay 
constants. 

Once quantum corrections are taken into account, loops of pseudoscalars can appear 
and the Lagrangian £4 starts to contribute. At next-to-leading order we have the three 
types of contributions depicted in Fig. [I], as well as wave-function renormalization. The 
diagrams are: i) the tree level diagrams giving a polynomial in quark masses and external 
momenta whose coefficients depend on the LECs of ( |2.2| ) via the diagram in Fig [l|(d). 
ii) the tadpole diagrams which contribute a logarithmic dependence on the quark masses 
from Fig. |l](b). In the present case these vanish because of the requirement ( |3.6| ) when 
combined with the contributions from wave-function renormalization. iii) the unitarity 
correction due to rescattering shown in Fig. |l|(c). 

We regulate using dimensional regularization and use the modified MS subtraction 
scheme as is standard in ChPT 0, Bl . 



Our results are in agreement with those presented in [12], [L3| and can be written in the 
form 



Fy(t)\ pi 



Ff(t) 



F 2 

TV 
1 



(2H(ml,ml,t) + H(m 2 K ,m 2 K ,t)) 



— (H(ml, ml,t) + 2 H(m 2 K , m 2 K , t)) 



F 2 

1 

F 2 



(-H(ml,ml,t) +H(m 2 K ,m 2 K ,t)) 



(3. 
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In terms of the integrals defined in App. [FJ, the function 7i is defined by 



1 2 

H(m 2 , m 2 , t) = -A(m 2 ) - B 22 {m 2 \ m 2 , t) + - tL 



(3.9) 



The form factors to 0(p 4 ) obviously satisfy the relation 



Ff(t) 



(3.10) 



At this order in the chiral expansion the expressions for the form factors depend on a 
single LEC, namely Lg, providing an opportunity for stringent experimental tests. This 
was discussed in detail in |12| with the then available data and will be discussed in Sect. 

B 



The choice of the pion decay constant F w in (|3.8| ) is a matter of choice. We could 
equally well have chosen F , the decay constant in the chiral limit, or F%. The difference 
is 0(p 6 ). 



3.3 Next-to-next-to- leading form factors 
3.3.1 General techniques and checks 

In order to obtain the expression for the form factors to 0(p 6 ) accuracy, one has to consider 
the diagrams of Figs. ^| and ^| in addition to wave-function renormalization contributions 
to the same order. At two-loops a whole set of new complications surface as compared to 
one-loop and especially the regularization and renormalization procedure has to be dealt 
more carefully. A long description tailored to ChPT can be found in Ref. |T3, and issues 
relevant to definitions of the subtraction are also treated in 



There are many ways in which to organize the calculation. One alternative is to use the 
one-loop expressions as generalized vertices and propagators in the two-loop calculation. 
This shortens the algebra but requires the knowledge of many one-loop processes also for 
off-shell external particles. We chose instead to perform explicitly the calculation of all 
Feynman diagrams in terms of lowest order quantities and afterwards put these into the 
physical quantities. A number of cancellations have to occur for both methods to agree 
and these, described below, form one of the checks on our calculation. 

Our calculations are made in dimensional regularization with the choice d = 4 — 2e. 
In order to include all finite contributions in d = 4, it is needed to keep loop integrals 
and constants in the diagrams of Fig. |2] to order e. In general any function, f(m 2 ,q 2 ), 
appearing in the loop diagrams has a Laurent expansion in e; 



f(m 2 ,q 2 ) = ^/ ( -^ m ^V-/ H) (™^ 2 )+/ ( ° ) K 1 ? 2 )+e/ (1) K ) 9 2 )+0(e 2 ). (3.11) 



;/' -'(//>7r,.'fH ' 

We never encounter terms more divergent than 1/e 2 and terms higher order in the expansion 
do not contribute to the order we work. The superscripts indicate the power of e the 
coefficient function gets multiplied with. As an example, diagram (d) in Fig. ^ contains 
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(a) 





(c) 











(f) 



(h) 
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(i) 



(j) 



00 



(1) 



Figure 2: The factorizable two-loop diagrams, one-loop diagrams with vertices from £4, 
and the tree level diagram with an £ 6 vertex, contributing to the form factors. The square 
stands for the insertion of the vector external current. A cross indicates a vertex from £ 4 
and a filled circle a vertex from £ 6 . 



two key ingredients: the one-loop two point function B (see Appendix ||) and the constants, 
L, appearing in £4, Eq. (f2.2|). Expanding both of them in Laurent series, one gets the 
following naive finite contributions 



diagram (d)| flnitc ai B^ l) L {1) + a 2 5 (0) L (0) + a 3 B {1) L { ~ 



1) 



(3.12) 



with <2j some polynomial that depends on the masses and/or external momenta. As shown 
in [20[, the constants can be chosen to be zero, their effect can be absorbed in the 
values of the 0(p 6 ) constants Cj. We will do so. The contribution containing B^ can be 
shown to always have counterparts cancelling it in other diagrams. This however requires a 
careful evaluation of all diagrams that makes this cancellation apparent and in particular, 
requires careful treatment of parts of the nonfactorizable diagrams. We have evaluated 
the latter in a more general fashion and we are thus obliged to keep the contributions 
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e 

(a) 



e 

(b) 



Figure 3: Non-factorizable two-loop diagrams contributing to the form factors. The square 
stands for the insertion of the vector external current. There exists also a crossed version 
of diagram (a), with the vector current insertion on the left hand side. 



containing 

Regularization and renormalization: 

In a local field theory, the final divergences are always polynomial in the masses and 
momenta. These can thus be subtracted from the freedom in the choice of parameters in 
the tree level Lagrangian. We use here the standard ChPT MS subtraction with e.g. for 
the coefficients of £4 

^ m "(Mh +i ' w ) (3 - 13) 

and 

c=-~(ln(47r)+r'(l) + l) . (3.14) 

Notice that Eq. (|3~T3| ) impl ies a choice of l/p as mentioned above. 

To 0(p 4 ) there are no problems and all divergences occurring diagram by diagram 
are polynomial. This is no longer true at 0(p 6 ) where there are two sources of nonpoly- 
nomial divergences. We have the one-loop diagrams with insertions of £4 vertices. The 
nonpolynomial part of the loop integrals is multiplied with ( |3.13 ) and leads to divergences 
proportional to ln(M 2 / /1 2 ) and others with momentum dependence. The required cancella- 
tion of these divergences leads to so called Weinberg relations. From the one-loop diagrams 
with £4 vertices we can determine the nonpolynomial divergences of the two-loop diagrams. 
These in turn determine the double poles. 

As an example of the above, the t dependent part of the nonlocal divergence produced 
by the diagram of Fig. |2|(f ) must be cancelled by other diagrams having the same possible 
divergence, in this case the diagrams of Fig. |(a,b,e,g,h,i) and of Fig. |3](b). The cancel- 
lation of both the t dependent and the ln(M 2 //x 2 ) dependent divergences is a very strong 
check on our calculation. 

Double logarithms: The cancellation as discussed above, can in fact be used to 
calculate part of the 0(p e ) correction using only one-loop diagrams. This has been done 



for quite a few processes in Ref. pl |. In the two- flavour case for many quantities these 



so called double-logarithm contributions often dominate, in particular they do so for mr- 



scattering and the corrections to the masses, |22|, [15| but not for e.g. the pion form factor 



[Hi. Our results here agree with the double-logarithm part as calculated in 
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Terms of the type Ll(ii) x L r j(fi): Similarly to the previous contribution the L^(ti) x 
L'j(fjL) pieces can be obtained in full generality without a complete calculation, since they 
only occur via tree diagrams and were also calculated in pi . Our results agree with those 
given there. 

Irreducible two-loop integrals: One of the main features of the non-factorizable 
integrals at two- loop (see Fig. fl3|)) in the three flavour case is that they can contain different 
scales even at zero momentum transfer. Those scales are given by the pseudoscalar masses. 
No exact analytic expressions are known for all the integrals we need here. Given the 
complexity of the integrals, it is worthwhile to have as many cross-checks on the numerical 
evaluation as possible. We differentiate between the two topologies in Fig. (|3]). The 
topology (a) is calculated with the methods described in 0. In particular, the subtraction 
constants are derived using the methods of [^3], E3, of which we obtained an alternative 
derivation of the recursion relations, and the dispersive part with a generalization of the 
methods of E5J to different masses. In addition we cross-checked the results using the 



techniques in |26], |27| . 

For the topology (b) we used the methods given in |2B], ^71] and cross-checked with \25\ 



for the case of equal masses. As an independent check for the tensor structures, we evaluate 
most of them numerically and then check that they satisfy the relations derived using the 
methods of Passarino and Veltman ||28|| , between the various different Lorentz structures. 

Putting the result in physical quantities: Up to this level we have computed all 
the diagrams in terms of the quantities appearing in the Lagrangian. In particular we have 
used the decay constant in the chiral limit and the lowest order meson masses. Once the 
amplitude is finite in terms of these bare quantities, with a cancellation of all nonlocal 
and nonpolynomial divergences as well as agreement of the polynomial divergences with 
the general evaluation of EO , we must rewrite these bare masses and decay constants 



in terms of the physical ones. As was mentioned already, the leading order in the form 
factors are mere constants fixed by Ward identities, Eq. (|3.6|) . We thus only need the 
0(p A ) expressions for masses and decay constants to do this rewriting. This is one of the 
reasons that the expressions here are so much simpler than in the case performed 
earlier, allowing us to quote the full result for the pion form factor, not only a numerical 
parametrization. 

There is of course a certain choice in how this writing in physical decay constants is 
done. At lowest order F w = F K = F v = F but at higher orders nothing prevents us to 
write the 0(p 4 ) expressions with Fjc or F v , rather than F n . 

For the masses the situation is easier in this case as well, the r] does not appear in the 
0(p A ) result, so the full 0{p A ) result can be simply written in terms of the pion and kaon 
mass. The masses can appear inside the arguments of the loop integrals. Rewriting those 
in the physical masses we expand the function in a Taylor series around the value of the 
mass at the next-to-leading order 

/Oio; <? 2 ) = firrif; q 2 ) + ^ (m| - m 2 ) Q^rf( m ^ OL^m, • (3.15) 
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For instance, for the scalar one and two-point functions, A and B, one obtains 



M™lo) = A{ml) + (ml - ml) B(ml, ml; 0) , 
B ( m lo, m mi Q 2 ) = B ( m l, m2 K> Q 2 ) + ( m lo ~ m l) C(ml, m 2 , m 2 K ; q 2 ) 

+ { m K0- m K) C(m 2 K ,m 2 K ,ml;q 2 ) , (3.16) 

where the m 2 are the bare or lowest order quantities and the m 2 are the equivalent renor- 
malized masses at next-to-leading order. Finally the C function is the scalar three point 
function 

dmlmlmU 2 ) = \ J d d P ^ _ ^ _ -L _ _ _ _ _ , ( 3 .17) 

Once these steps are performed, we are allowed to identify the masses and decay con- 
stants with the physical ones. Notice that this identification can be done without any 
manipulation in the genuine two-loop topologies from the very beginning because the dif- 
ference is of higher order. As a consistency check, we have cross-checked that all our final 
expressions are scale invariant. Notice that both the lowest order masses and decay con- 
stants and the physical ones are finite and this rewriting does not introduce any infinities. 

Three-point one-loop integrals: There is an extra check on the way that we per- 
formed the renormalization of the masses and decay constants. As is evident the diagrams 
(a,b,c,d,g) and (h) in Fig. (|2J) can be obtained by substituting the "bare propagator" of 
diagrams (b,c) in Fig. (JJ) by a "dressed" one-loop one. In this last method only one and 
two-propagator loop functions are permitted. In particular this tells us that the three- 
point functions, C, appearing naively via diagrams (g) and (h) should disappear once we 
promote all bare quantities to the physical ones. 

Gell-Mann— Okubo relation: There is one last comment concerning higher order 
corrections. At 0(p e ) we have chosen to rewrite all masses, except those inside logarithms 
and loop functions, in terms of the pion and kaon mass. The quark masses are written into 
those and the rj mass is removed using the GMO relation in the SU(2) isospin limit, 

3ml = Am 2 K - ml . (3.18) 

As mentioned above already, the vector form factors do not depend on the eta mass at 
next-to-leading order, only at next-to-next-to-leading order. This is not necessarily the 
case for other quantities, in particular the scalar form factors. Numerically at low transfer 
momentum the difference in using or not Eq. (|3.18 ) for the rj mass in the remaining loop 
functions is small. 

Independence of ji: A final check on the numerical programs is that they are inde- 
pendent of the scale \x when this scale is varied in the evaluation of the integrals and the 
values of the constants are changed appropriately. 

Comparison with previous work: As an independent cross-check of our calculations 
we compared whenever possible with the available expressions in the literature. This is the 



case for the K form factor which was obtained up to next-to-next-to- leading order in [17 
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for the isospin limiting case. We have compared partially our results for the same process 
and find agreement with the double logarithm terms and with the pieces containing the 



irreducible integrals. The same authors considered previously, [16], a certain combination 
for the vector form factor that leads to the Sirlin's relation f29fl . The q 2 dependence of 
this relation is free of unknown constants of 0(p 6 ) because of the Ademollo-Gatto theorem 
|30| . Some expressions for the pion form factor are also in the recent paper |jTj[] . The 
analytical comparison with this expression can only be made between the pieces obtained 
via the reducible diagrams finding agreement between them. For the rest of the expression, 



irreducible diagrams, the authors of [|16[ |18| use like us a numerical approach. The results 
of these are in good agreement. 

FORM: Most of the algebraic manipulations needed for this work were performed 
using FORM 1.0 and FORM 3.0 |31| while independent cross-checks of some structures 
were performed by hand. 

3.3.2 Results 

We split the full 0(p 6 ) contribution in several pieces 

F$(t)\ p6 = ij (F^(t) + F$ c (t) + Fy B (t) + F$ u (t) + Ffy(t)) , (3.19) 



FUt) 



with the following definitions: 
Fy L (t): The part which depends on the L\. 

The polynomial dependence on the Cq parameters C\. 
The part with the loop integrals that can be done analytically. 
The part that depends on the sunset integrals. 
The part containing the irreducible two-loop vertex integrals. 
The last three contain the proper two-loop diagrams and are the hardest part of the 
calculation. The precise separation between them is somewhat dependent on how the set 
of nonfactorizable two-loop integrals is chosen. In addition, there are rather large numerical 
cancellations between these three parts because the method to separate factorizable from 
nonfactorizable contributions often induces large terms in both which afterwards cancel. 

We have not performed the separation into tadpole and unitarity corrections here. Be- 
cause of the requirement ( |3.6| ) the tadpoles cancel to a very large extent and the separation 
between tadpole and unitarity contributions is in any case representation dependent. 
The dependence on the L\ is fairly simple: 

Fv L (q 2 ) = %ml{2Ll + Ll)A{ml)+AmlLlA{ml) + q 2 Ll{QA^ 

+B 22 (ml, ml; q 2 ) {-16m 2 {2L\ + L\) + 4g 2 (4L^ - 2U 2 + 2U 3 - U 9 )j 
+B 22 (m 2 K ,m 2 K ;q 2 ){-8mlLl + 4q 2 Ll-2q 2 Ll} , (3.20) 



^vlV) = (l6m^ + 8m^)A(m^)+4m^(m 2 ) + g 2 ^(5A(m 2 )+4A(m^)) 
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Figure 4: The real and imaginary parts of the loop diagrams at 0(p 4 ) and 0(p 6 ) with all 
L\ = 0, for the pion form factor. Notice that there is convergence both for the real and 
the imaginary part. 

+B 22 (m 2 K , m 2 K - q 2 ) {-32m 2 K L r 4 - 16w^ + Aq 2 {AL\ - 2L r 2 + 2L r 3 - U)} 
+B 22 {ml ml- q 2 ) {-SmjLg + - 2q 2 L r 9 } 

16U 5 Llq 2 (ml - m 2 K ) , (3.21) 

and 

F$l{q 2 ) = {AmlLl + q 2 Ll) {-A(ml)+A(m 2 K )) 

+B 22 (ml, ml; q 2 ) {8m 2 ^ - 2q 2 {2Ll - U g )j 

+B 22 (m 2 K ,m 2 K -,q 2 ) {-8mlLl + 2q 2 (2Ll - L r 9 )} . (3.22) 

We can now check how large the various loop contributions are. In Fig. [| we have 
shown the real and imaginary parts of the pure loop contributions at 0(p A ) and 0(p 6 ) for 
the pion form factor. To put the figures in perspective one should bear in mind that the 
lowest order is exactly one and that at t — ±0.25 GeV 2 the 0(p 4 ) contribution from Lg 
setting Lg = 6.9 x 10~ 3 is ±0.40. We show the difference between the contribution from 
the pure loops only for the three different pseudoscalars in Fig. |5] at 0(p A ) and in Fig. |6| 
at 0(p 6 ). Fig. shows the contribution from the L\ dependent part at C(p 6 ), Fy^it) 
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loops only l_[ = C[ = 



u_ 




-0.25 -0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 0.25 

t [GeV 2 ] 



Figure 5: The real parts of the loop diagrams at 0(p 4 ) for the pion, charged kaon and 
neutral kaon form factor. 



with the input parameters as given in pH 



The part that depends on the Cq constants C\ is also rather short. 

Fv C (q 2 ) = -4q 2 ml(4C r 12 + 4Cl 3 + 2C; 3 + C r m + C r 65 + 2C; ) 

-8q 2 m 2 K (4C{ 3 + C r M ) - Aq\C; 8 - C r 90 ) , (3.23) 



*veV) = -q 2 < (l6C[ 3 + |cS, + 4CT M - ^ 

- q V K (lQC{ 2 + 32C[ 3 + y C r 63 + 8C* M + y CJb + 8CSo 
-4g 4 (C^ - CJ,) , (3.24) 



and 



^c(g 2 ) = 3? 2 K-0(2Q 3 -C65)- (3-25) 
At low-energies, for instance the space-like interval with y/—t < 350 MeV, the pion 
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loops only l_[ = C[ = 

0.1 
0.08 
0.06 



0.04 

Q-> 

0.02 


-0.02 
-0.04 

-0.25 -0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 0.25 

t [GeV 2 ] 

Figure 6: The real parts of the loop diagrams at 0(p 6 ) with all L\ = 0, for the pion, 
charged kaon and neutral kaon form factor. The scale is the same as in Fig. [|. 

vector form factor can be expanded in a Taylor series as 

F v {t) = 1 + l -(r 2 Y v t + c v t 2 + 0(t 3 ) , (3.26) 
o 

defining in that manner the radii, (r 2 ) v , and the coefficient c v . 
3.4 Isospin breaking 

Throughout this paper we work in the isospin limit. The isospin breaking contributions 
can be split naively into two parts. 

A first contribution comes through electromagnetic interactions. This has been studied 
for the pion form factor in f32fl . This is primarily obtained by taking photon loops that 
develop an infra-red singularity. Near threshold in the time-like region this contribution 
can dominate over the next-to-next-to-leading terms because of Coulomb final state inter- 
actions. This is not the case for the pion form factor itself, the corrections are about 1%. 
For the constant c v there might be a large contribution coming from resolution dependent 
terms but they cancel to a large extent for typical values of the detector resolution. The 
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Figure 7: The real parts at 0(p 6 ) for the pion, charged kaon and neutral kaon form factor 
that depend on L\. This corresponds to Fy h (t). The inputs are taken from |I0[ . The scale 
is the same as in Fig. |5|. 

electromagnetic contribution to the masses we take into account indirectly, we use the 
charged pion and charged kaon mass in our numerical analysis. 

The second source of isospin breaking is from the mass difference m u — ma- The pion 
form factor receives corrections starting only at O {{rrid — m u ) 2 ) [RD|. This together with 



the constraint ( |3.6|) means that corrections only come from loop diagrams and are thus 
expected to be small. In the kaon form factors isospin breaking could be larger but given 
the present experimental accuracy is not expected to be relevant. 



4 Data and fits with the ChPT expression 
4.1 Data on Fy 

Data on the pion form factor have been collected in both time-like, t > 0, and space-like, 
t < 0, regions. We discuss both in turn. 
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Data 


n ( r2 )v [f m2 




FNAL H 
NA7 1 34 j 


0.973 ± 0.033 0.382 ± 0.078 
= 1 0.443 ±0.018 

0.990 ±0.004 0.428 ±0.011 
= 1 0.452 ± 0.006 



Table 1: The parameters of a simple pole fit to the ire scattering data. Note the difference 
with the normalization free. 



4.1.1 Space-like 

In the space-like region the measurements have been done using ire — > ire scattering. There 
are two main experiments [ 33 1 at FNAL and NA7 at CERN p4 |. Both experiments are 



based on pions scattering off the electrons of a liquid hydrogen target and the data agree 
in the overlap region. The last one has significantly smaller error bars and dominates 
clearly all the fits. There are also somewhat older FNAL data from the same experiment 



using a lower energy pion beam [35]. These have larger errors and are compatible with the 
newer data. References to older data can be found in both [ 34" ] and J33J . These data were 
analyzed in both papers with a pole-fit 



n 



2 ' 



where n refers to the normalization uncertainty. Ref. 



obtained 



(r 2 )^ = (0.431 ± 0.010) fm 2 

when they constrained n and 

(r 2 )^ = (0.427 ±0.010) fm 2 
with n left free. The equivalent results of [[33 



when they constrained n and 



are 

(r 2 )^ = (0.439 ± 0.030) fm 2 



r 2 y y = (0.384 ± 0.088) fm 2 



(4.1) 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



with n left free. The two data sets together with the results of our own pole fits are shown 
in Fig. H and the pole fit parameters are in Table [l]. The difference with the numbers 
quoted by the experiments themselves is because there the normalization was fixed using 
the experimental error rather than set exactly equal to one. The quoted errors for the 
normalization are 1% for [33] and 0.9% for [|54j|. 

Using the previous data sets there are several determinations of the pion vector form 
factor in the literature. The experiments were analyzed directly mainly with a pole fit as 
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7te scattering data and pole fits 




0.4 - 



0.3 ^ 1 1 1 1 1 

-0.25 -0.2 -0.15 -0.1 -0.05 

t [GeV 2 ] 

Figure 8: The two data sets of ne scattering and the resulting pion form factor squared. 
Also shown are pole fits to both data sets. 



discussed above. A full analysis within the framework of two flavour ChPT is in Ref . |L5| , 
and we will compare with the results from there in Sect. |5|. 

In [[34] the data were also analyzed using a dispersion relation technique with a Pade 
approximation for the Omnes representation and irn data for the 8\ phase shift the result 
they obtained is 

(r 2 )^ = (0.439 ± 0.008) fm 2 . (4.6) 
Details of this technique can be found in [j36T]. 



4.1.2 Time-like 



In the time-like region there are several experiments, but unfortunately most of the data is 
located outside the range of applicability of ChPT. Essentially the data is obtained from 

j, 4lJ . In the later experiment, [fiT 



-R-RU-r 



37 and e e H 



7T + 7T 



38, B9 



two sources: r 

corresponding to the NA7 collaboration, the systematic errors are close to the statistical 
ones and we believe that they have been underestimated in the time-like region as they are 
difficult to fit together with the space-like data. 

A resonance based model was used in to fit the data up to 1.2 GeV and determines 
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the mass and decay constants of the p as well as the pion radius with the result 

(r% = (0.460 ± 0.011) fm 2 . (4.7) 
Instead one can rely on r — > ixixv T data to fit the expression 

dz 5\ 



f 1 t 3 f°° 

F*(t) = exp < ait + -a 2 t 2 + - / 

{ 2 71 J im 2 



z 3 z — t — it 



(4- 



with ai,a2 given in terms of (r 2 ) v and dy (see |42|]) and 5\ determined by resonance 
saturation. Keeping only the p contribution one obtains 



r y v = (0.429 ±0.01) fm 2 . (4.9) 



4.2 Data on F v 



K+ 



The charged kaon form factor has been measured directly via kaon electron scattering 
method by the same two main experiments that measured the pion form factor. The 
older experiment is a FNAL experiment JEJ and reported a rather low kaon charge radius 
obtained by a fit to a pole form factor 

(r 2 )£ + = (0.28 ± 0.05) fm 2 . (4.10) 

They also quoted a dipole fit where the normalization was left free with the result 

(r 2 )£ + = (0.42 ± 0.20) fm 2 . (4.11) 



In [|33| a simultaneous analysis of the pion and kaon form factors was performed from the 



obtained data quoting the difference 

(r 2 }y + ~ (r 2 ) v = (-0-16 ± 0.06) fm 2 . (4.12) 

This difference is claimed to be largely free of systematic errors. 
The more recent CERN experiment |44j] finds using pole-like fits 



(r 2 )£ + = (0.34 ± 0.05) fm 2 , (4.13) 
and from the ratio to the simultaneously collected pion data 

(r 2 }y + ~ (r 2 )v = (-0-10 ± 0.045) fm 2 . (4.14) 

Both data sets are shown in Fig. || together with a simple pole fit with the normalization 
fixed to 1 and left free to both data sets. The parameters of those fits are in Table 0. 

The threshold for kaon production in tau decays or e + e~ annihilation is too high for 
the ChPT expressions to be valid. In principle a dispersion relation could be used to relate 
those data with the threshold parameters but this involves a rather large extrapolation. 
We have therefore not included this type of data in our analysis. 
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Ke scattering data and pole fits 




-0.12 -0.1 -0.08 -0.06 -0.04 -0.02 

t [GeV 2 ] 

Figure 9: The two data sets of Ke scattering and the resulting kaon form factor squared. 
Also shown are pole fits to both data sets. 



4.3 Data on Ff 

The data on the neutral kaon electromagnetic form factor are rather scarce. There is a 
measurement of the radius using kaon regeneration on electrons [[EJ leading to 



(r 2 )£ = (-0.054 ± 0.026) fm 2 . (4.15) 

The decay Kl — > 7r + 7r~e + e~ also has contributions that contain the neutral kaon electro- 
magnetic form factor but the extraction from the data is rather model dependent and not 
yet done at present. 



4.4 Input parameters 



For the physical masses we use the charged pion and kaon mass |46| and the pion decay 
constant, 

m w + = 139.56995 MeV , m K + = 493.677 MeV and F n = 92 A MeV . (4.16) 

For the eta mass we use the mass resulting from the GMO relation. As discussed above, 
this does not make any difference within the accuracy of our results. 
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Data 


n ( r2 )v + P™ 2 




FNAL H 
NA7 44] 


1.084 ±0.091 0.457 ±0.198 
= 1 0.278 ± 0.046 

1.019 ±0.030 0.395 ±0.109 
= 1 0.334 ± 0.045 



Table 2: The parameters of a simple pole fit to the Ke scattering data. Note the difference 
with the normalization free. 



A somewhat more complicated issue is how to deal with the input values of the L[ other 



than Lg. We use the two-loop accuracy determination given in fit 10 of [10]. We choose fit 
10 and not the main fit, because the analysis of the E865 experiment has been confirmed 
in the meantime ||47||P]. Notice that in doing so, i.e. using the values in |10[ without any 
new refitting, we have assumed that the value of Lg does not have a substantial impact on 



the rest of LECs. This is true, as we discussed in ||. For the quantities studied in |L0 
Lg only appears inside the form factors, but together with the squared effective mass 
of the dilepton system and varying Lg over an extremely wide range did not change the 
fitting results. 

4.5 ChPT fits to F$(t) 

Checking the Lg dependence of Fy, we find that the real part of this can be well described 
by a polynomial of the form 

Fv(%e >L r 9 =L r 9 {at + bt 2 + ct 3 ) . (4.17) 

The coefficients a, b and c are in principle dependent on the other L\ used as input. 

We fit the data on Fy(t) with the ChPT expressions with the value of Lg set equal to 
zero and a polynomial of the form 

a f t + b f t 2 + c f t 3 (4.18) 

for several variations of the input. In this way we can separate the dependence of Lg on 
the unknown 0(p 6 ) constants. We fit similarly at 0(p 4 ). 

In the fit we have left the normalization of the data from [33], [Q and [51] as an 
additional fit parameter. 

The fits we have performed are a pure 0(p 4 ) fit, an 0(p 4 ) fit but allowing also a t 2 and 
t 3 polynomial, and the same for the 0(p 6 ) expressions. For the input of the L\ we used 
fit 10 of ||10|| . The fits to the other L\ with varying the assumptions on L\ and Lg from 
that reference make essentially no change as can be seen from the rows labelled fit 11, fit 

1 Note that we only use the linear fit of the experimental slopes and derivatives at threshold for the form 
factors and not their quadratic fit whose curvature can not be reproduced by the next-to-next-to-leading 
calculation in ChPT of the Kn form factors with reasonable choices of the parameters. 
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n 



0( P 6 ) fit 11 
0(p 6 ) fit 12 
0(p 6 ) fit 13 
polynomial only 

0(p 6 ) space 



3| npf nfl 



a f [GeV]- 2 6, [GeV]~ 4 c f [GeVf 



0.981 
0.995 
1.001 
0.994 
0.994 
0.994 
0.994 
0.995 
1.000 



0.979 
0.998 
1.009 
0.998 
0.998 
0.998 
0.998 
0.998 
1.005 



1.073 
0.862 
0.871 
0.860 
0.860 
0.860 
0.860 
0.871 



1.285 ±0.003 
1.647 ±0.054 
1.788 ±0.045 
1.693 ±0.054 
1.696 ±0.054 
1.704 ±0.055 
1.690 ±0.054 
1.936 ±0.056 
1.844 ±0.018 



= 
3.83 ±0.28 
3.16 ±0.27 
3.45 ±0.28 

3.44 ±0.28 
3.43 ±0.28 

3.45 ±0.28 
4.34 ±0.28 
5.53 ±0.16 



= 
5.98 ± 1.69 
= 

5.33 ± 1.69 
5.32 ± 1.69 
5.30 ± 1.72 

5.34 ± 1.69 
6.44 ± 1.69 
12.3 ± 1.0 



Table 3: The results of fitting to Fy(t) with different theoretical inputs as explained in the 
text. = a stands for a value held fixed equal to a. 

12 and fit 13. We also presented a pure polynomial fit. The x 2 for all the fits including 
the t 3 term is the same and the quality of the fit is good. The x 2 is 60.3 for 75 degrees of 
freedom. The last line indicates the impact of the time-like data. This fit was performed 
with the space-like data only and had a x 2 = 48.6 for 54 degrees of freedom. 

Notice that the fits are somewhat different, but compatible within errors, with the ones 



given in The reason for this difference is that the fits there were performed with the 



data of 41] with the normalization fixed to one. We have left that free. 



4.6 ChPT fits to Ff + 

We now perform the same fits to the kaon data. The time-like data are very much outside 
the domain of validity of ChPT so they are not included. The space-like data are much 
less precise than the pion ones and have a much smaller range of t. There will thus be 
much less information about the higher powers of t. For the pion we could fit using the 
whole region —0.25 < t < 0.25, while here we have only the range —0.11 < t < —0.02. 
More care is needed in performing the fits, for instance, fitting with a t 3 term and both 
normalizations free leads to rather nonsensical results as can be seen in the sixth row of 
Table |. 

4.7 Results for the charge radii and Cy 

The NLO ChPT expression for the pion form factor [|l|, [13| was already analyzed quite 



some time ago, using the space-like data ||13fl . The results were 

(r 2 ) v = 0.392 (0.366) fm 2 , (4.19) 

with a fit with normalization one (free). 

The analysis using two flavour ChPT at two-loops led to the conclusions fTS] 

(r 2 ) v = (0.437 ± 0.016) fm 2 , c\ = (3.85 ± 0.60) GeV 2 (4.20) 
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Ti 


441 


77 


43; 




Of \GeV]~ 2 


bt fGeVl -4 


ct [GeVl~ 6 

o/ [ore v j 




I. 


015 


l. 


055 


1.408 ±0.300 


= 


= 




0. 


993 


l. 


02£ 


\ 


0.897 ± 1.169 


-5.80 ± 11.60 


= 


0(p 6 ) 


1. 


015 


l. 


055 


1.464 ±0.299 


= 


= 


0(p 6 ) 


1. 


027 


l. 


070 


1.756 ±0.297 


= 3 


= 


polynomial only 


0. 


890 


0. 


919 


-3.259 ±3.989 


-96.9 ±82.0 


-559 ± 492 


polynomial only 


0. 


993 


l. 


02£ 


\ 


0.978 ± 1.169 


-5.51 ± 11.60 


= 


polynomial only 


1. 


014 


l. 


054 


1.515 ±0.300 


= 


= 



Table 4: The results of fitting to Fy + {t) with different theoretical inputs as explained in 
the text. 



including the data of [41| and 

(r 2 y v = (0.453 ± 0.016) fm 2 , = (4.45 ± 0.60) GeV 2 

without including them. 

The polynomial fit to all data gives 

(r 2 )l = (0.452 ± 0.013) fm 2 , c\ = (4.34 ± 0.28) GeV 2 

The pure 0(p 4 ) fit yielded 

(r 2 )^ = (0.368 ± 0.001) fm 2 , 
while the 0(p 6 ) fit including the t 3 term leads to 

(r 2 y v = (0.452 ± 0.013) fm 2 , c\ = (4.49 ± 0.28) GeV 2 



(4.21) 



(4.22) 



(4.23) 



(4.24) 



This last number we consider our final result for these quantities. We have reproduced the 
previous fits with similar underlying assumptions. 

The kaon charge radius obtained from the linear fit is 



r 2 )* + = (0.354 ± 0.071) fm" 



(4.25) 



while fitting the 0(p A ) expression leads to 



(r% = (0.361 ±0.071) fm 2 



(4.26) 
(4.27) 



and the 0(p 6 ) fit to 

(r 2 )£ + = (0.363 ± 0.072) fm 2 . 

Finally, fixing the quartic slope to a value similar to the one obtained from the pion, we 
obtain 

(r 2 )£ + = (0.431 ± 0.071) fm 2 . (4.28) 



23 



5 Discussion and Determination of Lg 

The ChPT expression for the form factors depends on the unknown constants CJ". For the 
remainder of the discussion we define 

= <<z 2 + < 2 Y, (5.1) 

and the only model-independent relations are 

Rv2 = Rv2 and Rf 2 = . (5.2) 

5.1 Estimates of the Ry i 

In this section we estimate some 0(p 6 ) constants. In principle these constants should be 
obtained from QCD directly but this is not possible at present. This, together with the 
large number of LECs at higher orders in the chiral expansion, makes it unavoidable to 
resort to model estimates to gain some predictability. We stress that this model dependence 
only starts at 0(p 6 ). At 0(p 4 ) we use the data directly. The main idea is to saturate the 
0(p 6 ) order parameters via the exchange of higher mass resonances HR fEJ). We use the 



notation from |48j to which we refer for a more extensive exposition of the method. One 



uses a matter field Lagrangian coupled to the Goldstone Bosons octet 

£ rcs = Yl { £ k,„(^) + &«} , (5-3) 

V,A,S,P 

where for the cases we are interested in, the allowed intermediate states are reduced to 
vectors. For these spin-1 mesons we use the realization where the vector contribution to 
the chiral Lagrangian starts at 0(p 6 ). We shall specifically only discuss terms relevant for 
our purposes. 

/ ' /y :(^K ) ^]) + / x (^K, X _]}, (5.4) 



2^2 



with 



and Vfj, is a three-by-three matrix describing the full vector nonet. Furthermore my refers 
to the octet and singlet masses in the chiral limit. The matrix is defined in terms of 
the covariant derivative, Eq. (|2.5|), 



= iw D^Uw = vl , u 2 = U . (5.5) 
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Instead of evaluating each contribution diagrammatically we integrate out the heavy 
degree of freedom by means of the classical equation of motion obtaining 



-|^<(v A ^0(v> M ,M)> 



fxfv 



<(v A /f)K,x-]> 



y/2>, 



mi 



(5.6) 



where the resulting Lagrangian is directly of 0(p 6 ). 

The numerical values for the couplings constants appearing in ( |5.6|) are obtained using 



experimental information. The value of fy is obtained from p — > e + e |4"5fl. The parameters 



f x and gv are determined simultaneously from p 



1 A 

we use 



7T7T 



and K* — » fT7r ]T9 ] . In summary 



/v = 0.20, A 



-0.025, V = O.O9. 



and the vector- mass is taken as the experimental one [4(| 



my = m. 



0.77 GeV. 



(5.7) 



(5.* 



The description is basically identical to the one in |15| but now also includes the kaon case. 
Numerically we obtain 



D7T 

R Vl 

R K+ 



R 



VI 

K° 
VI 



= 0.26 x 10~ 3 , 
-0.49 x 10~ 5 GeV 2 
-6 x 10~ 5 GeV 2 , 
0. 



(5.9) 



An alternative estimate is based on the full VMD form for the form factors 

,r 2 

Fyit) 



m r 



t 



FT® 



Ff(t) 



m r 



2 m 2 



+ 



m 



+ 



t 6 m 2 , — t 3 m 2 — t 



1 m l . 



+ 



1 m 



+ 



2 m 2 —t 6 m 2 , —t 3 m 2 , — t 

p UJ w 



With m v as above the chiral limit of the vector meson masses, this leads to 



(5.10) 



T}TT 

tiv 



V2 



K V1 



r>K+ 
K V1 



tiy\ 



r>K+ 

n V2 



0.21 x 10- 



rn. 



v 



m 



m 



v ) « -0.35 x 10~ 5 GeV 2 



v 



m 



v 



ml 



1 2 , 1 2 , 1 2 2 

2 P 6 w 3 ^ y 



1 2 i 1 2 i 1 2 2 

—m n H — m, , H — m, — mv 
2^63 



-3 x 10~ 5 GeV 2 , 
-3 x 10~ 5 GeV 2 



(5.11) 
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We used here m? p — rriy « (m 2 K * — mV)rri^/rn 2 K . The pion numbers of ( p . 1 1|) and (O) 



are in reasonable agreement. The numbers for the kaons have a sizable uncertainty. E.g. 
changing ray to in the denominator changes the values by a factor of (m p /m^) 4 « 0.3, 
which differs from one only by higher order effects. We will use the values for Ry 1 and 
R V2 but not for the kaons, sizable higher order quark mass effects are known to exist in 



the vector meson sector 50 



5.2 Obtaining L\ 

At 0{p A ) the relation between a/, Eq. (|4.18|) , and L\ is direct, 1L\jF 2 = aj and we obtain 

L£(0.77 GeV)| p4fit = 5.5 (7.0) x 10" 3 . (5.12) 

The number in brackets is from the fit with bf and c/ left free. The value obtained from 
the kaon form factor is perfectly compatible with this since the values of a/ are compatible 
within errors. 

At the next-to-next-to-leading order the relation is a little different 

a f = ^ (1 + Ag) L r 1 ^ (513) 

7T 7T 

The value of A 9 follows from our calculation and depends explicitly on the values of the 
other L\. For the inputs of fit 10 in |H| we obtain 

A 9 = 0.27. (5.14) 

This together with the value of Ry l of ( |5.9| ) leads from the 0(p 6 ) fitted value of a/ to 

LgXO.77 GeV)| p6fit = 6.25 (5.93) x 10~ 3 , (5.15) 

where the number in brackets is the fit including the t 3 term. The difference between the 
two numbers is an indication of the 0(p 8 ) and higher corrections. Neglecting R vl would 
have lowered both these numbers by 



L^O.77 GeV)|^ = 0.23 x 10^ . (5.16) 



Taking the difference in the two values of L\ in ( |5.15| ) as an estimate of the error from 



higher orders, allowing for a factor of two uncertainty in R vl and adding the experimental 
error, all in quadrature leads to 

L r g (0.77 GeV) = (5.93 ± 0.43) x 10~ 3 (5.17) 

as our final result. Notice that the error due to experiment only is about half this. 

Using this value for L r 9 and the value of bf we can extract a measured value of R V2 of 

R* 2 = -4 (C 8 r 8 - C; o ) = (0.22 ± 0.02) x 10" 3 , (5.18) 



in good agreement with the estimates ( |5.9j ) and ( |5.11| ) 



The kaon data do not lead to any better determination of these parameters. The extra 
constraints we obtain from them are given in Section |5T3", 
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5.3 Charge radii and the Cy in ChPT 

We now present the different contributions to the charge radii. For the pion charge radius 
the different contributions are given by 

(r 2 ) v = |o.325[L 9 p 4 } + 0.068[loops p 4 ] + 0.084[L 9 p 6 ] 

-0.01 1 [loops p 6 ] -0.015^]} fm 2 
= 0.451 fm 2 , (5.19) 

where the terms inside square brackets reveal the source. The charged kaon charge radius 
gives similarly 

(r 2 )£ + = |o.325[L 9 p 4 } + 0.031 [loops p 4 ] + 0.018[L 9 p 6 ] 

-0.011 [loops p 6 } - 0[Rvt]} f m 2 
= (0.363 ± 0.19) fm 2 . (5.20) 



The error is given by using the value of ( |5.9| ) for Ry x . We can in fact use these results 
to obtain an experimental bound on the combination of 0(p 6 ) constants containing terms 
with masses 

(r 2 )v ~ (r 2 )f = 0.102 fm 2 + A (ifo - <") . (5.21) 

7T 

With the experimental value (|4.14j) we obtain 

R V1 - R^ = (0.0 ± 1.4) x 10" 5 GeV 2 (5.22) 



or with the experimental value Q4.12 ) 

R^ - R^ = (1.9 ± 1.9) x 10~ 5 GeV 2 . (5.23) 



Both values are within a factor of two to three of the estimates done above in (|5.9|) and 

(pm 



Our calculation can be used to predict the neutral kaon charge radius. The Lg depen- 
dence is basically zero at 0(p 6 ) as well. We thus obtain 

(r 2 )f = | - 0.0365 [loops p 4 ] - 0.0057[loops p 6 ] } fm 2 + ^Ryl 

7T 

= (-0.042 ± 0.012) fm 2 . (5.24) 

The error is based on assuming that the unknown contribution is not larger than twice the 
0{p & ) loop contribution. Using the estimate for Ryl of ( |5.11| ) lowers the value of ( |5.24| ) 
by about 0.10. The result is in good agreement with the measurement ( f4.15|) . We can also 
turn the argument around and obtain 



R V2 = (-0.4 ± 0.9) x 10~ 5 GeV 2 , (5.25) 
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which is compatible with the estimate ( |5.9| ) and somewhat lower than ( |5.11[ ). 



The contribution to c v can be expanded similarly. For the pion we obtain 

cl = |o.66[loops p 4 } + 0.41 [L 9 p 6 } 

+0.38 [loops p 6 } + 3.04[i^ 2 ]j GeV" 4 

= 4.49 GeV~ 4 . (5.26) 

Here we used the estimate for R V2 obtained from the measurement (fifth row bf — 3.45) 
in Table ^| and removing the contribution from Lg. For the charged kaon 

4 + = |o.37[loops p 4 } + 0.19[L 9 p 6 ] 

+0.11 [loops p 6 } + 3.04[#{£]} GeV" 4 
= 3.71 GeV^ 4 . (5.27) 

Note that the data as discussed are not good enough to give a reasonable value for this 
parameter. Finally we obtain for the neutral kaon 

cf = | - 0.29[loops p 4 } - 0.22[L 9 p 6 } 

-0.23[loops p 6 } + 0.00[i^ 2 + ]} GeV~ 4 
= -0.74 GeV~ 4 . (5.28) 

We also show the overall agreement with the data in Figs. 10 and for the measured 
pion form factor. Note that we have plotted the data with normalization one, not the fitted 
normalization. As can be seen the convergence is nice and of a similar quality as the two 
flavour results in |~5|. As stated before, we are in excellent agreement with that reference 



when the difference in the treatment of the data of j|T| is taken into account. 



5.4 Comparison with predictions for Lq 

Our estimates for the 0(p 6 ) are based on a resonance saturation model. As seen above, 
unless the estimates are off by more than a factor of two, the value obtained for L r 9 is 
reliable. In table we compare the result fl5.15| ) with several model approaches. First of 



all with the resonance saturation prediction |48| that in this particular case reduces to the 
Vector contribution 

Ll{mv) = ^f- (5 - 29) 

In addition if one implements some extra "QCD inspired assumptions" for the high energy 
behaviour, such as an unsubtracted dispersion relation for the pion form factor, one obtains 

L5(m v ) = ^-. (5.30) 
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full result 




0.3 - 

0.2 ' — 1 1 1 1 1 

-0.25 -0.2 -0.15 -0.1 -0.05 

t [GeV 2 ] 



Figure 10: The comparison of the space-like measurements with the ChPT calculation. 
Notice that there is excellent convergence over the whole kinematical range. 



We also show a constituent-quark model |51J which in the chiral and large-N c limit and 
including leading gluon contributions leads to 



48vT 



:[1 + 0(1/M$ 



(5.31) 



with Mq being a constituent quark mass. Using only the first resonance plus a continuum 
spectrum to mimic QCD one obtains the so call LMD model [ 5"2"|| . If one implements in 



addition QCD sum rules, one gets f53 



[m. 



in 

2 m 2 p 



(5.32) 



An extensive discussion can be found in the more recent study of ||54 |. At the expected 
accuracy all of these predictions are in good agreement with the value obtained in ( |5.17| ) . 



6 Summary 

Let us summarize our findings. We have described the evaluation at 0(p 6 ) in ChPT of the 
pion and kaon electromagnetic form factor, giving an overview of the methods and checks, 
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full result 




0.05 0.1 0.15 0.2 0.25 



t [GeV 2 ] 



Figure 11: The comparison of the time-like measurements with the ChPT calculation. 
Notice that there is excellent convergence over the whole kinematical range and the slow 
buildup of the tail of the p-meson. The data have not been corrected for the normalization. 

both numerical and analytical. The pion vector form factor expression is quite handleable 
and we have listed it completely, partly in App. fA] and the rest in Sect. |3.3.2 . 



We have assumed a resonance saturation model to determine the Born contribution 
appearing at 0(p e ) and compared it with the 0(p 6 ) parameters we could determine directly. 
There is a priori no reason to doubt that this or similar methods retains the bulk of the 
0(p 6 ) counter terms. This point was certainly fulfilled for all the studies inside the two- 
flavour case. What is evident is that the use of an effective approach enforces to use an 
estimate for the higher order parameters. This constitutes at present the main source of 
theoretical uncertainty. 

If we start with the premise of a "well behaved" series expansion such terms are sub- 
sub-leading and therefore the impact of their precise value ought to be mild. Thus, to make 
some sense of the full approach, there should be no need to fine tune many of the 0(p 6 ) 
constants. In this work, the impact of these higher order parameters seemed reasonable. 

For the pion form factor we have collected all the available data with reasonable pre- 
cision and combined it with our analytical results to obtain a new determination of the 
pion charge radius and dy, as well as the LEC Lg(m p ), Eq. (|5.15|) . All these results are 
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Eq. (|5.15) 


48 


|52|, 




51| 


5.93 ~ 7 6.5 6.33 



Table 5: Values of 10 3 x Lg(m p ) for different model approaches in comparison with its 
value obtained with an 0{p & ) fit. 

compatible with previous ones. In addition we have evaluated the correction due to higher 
order terms. The genuine loop contributions do not show any strong deviation and they 
define a proper convergent expansion, see Fig. [|. 

The data on the kaon form factor are also included and are within errors perfectly 
well described by ChPT. We have used these data to put experimental bounds on two 
combinations of 0(p 6 ) constants involving quark masses. 

It is also evident from our discussion that the direct role of the Zweig rule suppressed 
terms with L\ are marginal for the present processes because they are accompanied by 
products of m 2 . The indirect impact via its influence on the determination of the other L\ 
fT(| is also small here. 
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A Pion vector form factor 

In this appendix we give the explicit expressions for the 0(p 6 ) parts of the pion form factor 
which were not given in Sect. |3.3.2| . 

A.l Reducible contributions: 

1-loop contribution x 1-loop contribution 

This contribution can be obtained from diagrams (a,c,e,g,i) and (j) in Fig. In terms of 
the integrals defined in App. [B] we obtain 

Fvb(Q 2 ) = l/(l6n 2 ){-l/2mlA(ml)-l/4mlA( m 2 K )} 

+1/(16tt 2 ) 2 {-1/48 ml{3m 2 K + lOm 2 )^ 2 - 35/96 m^(2m| - m 2 ) 
-89/48 ml - l/16g 2 (m^ + 2m 2 )(l + tt 2 /6)} 
+ l/(167r 2 )^ 2 (m 2 , ml , q 2 ) (5m 2 - 1/2 q 2 ) 
+ ini6n 2 )Bl 2 (m 2 K ,m 2 K ,f)(3/2ml-l/^q 2 ) 
+4B 22 {m 2 n , ml, q 2 ) 2 + 4 J B 22 (m 2 , m 2 , q 2 )B 22 {m 2 K , m 2 K , q 2 ) 
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+B 22 (ml,ml,q 2 ) {-4A(m 2 ) - 2A{m 2 K )} + B 22 {m 2 K , m\, q 2 ) 2 

+B 22 {m 2 K ,m 2 K ,q 2 ) {-2A(mD -A(m 2 K )} 

-1/4 A(m 2 ) 2 - 3/8m 2 Jm 2 K A(m 2 K ) 2 + A(ml)A(m 2 K ) 

+ l/4A(m 2 K ) 2 - 1/8 q 2 /ml~A{ml) 2 - 1/16 q 2 / m 2 K A{m 2 K ) 2 . (A.l) 

A. 2 Irreducible contributions 

We have disentangled the contribution to the irreducible diagrams in two pieces. The first 
one comes from the topology (a) shown in Fig. |3|. The treatment of the possible tensor 
integrals was treated lengthly in J7| where we refer the reader for the notation. This term 
reads 

*Sn= H F '(mlmlml;ml)(3/2mt) 

+H F '(ml,m 2 v ,m 2 v ,ml)(l/18mi) 

, ttF' I 2 2 2 2\i 2 2 \ 

+H [m K ,m^,m K ;m m )(m^m K ) 

+H F '(m 2 K , m 2 K , m 2 ; m 2 )(-5/6 m*) 

+H F \m 2 ,m 2 K ,m 2 K ;ml)(l/Aml(2m 2 K -l/2ml)) 

+H['(ml, ml, ml; m 2 )(-2rr4) 
+Hr(ml,m 2 K ,m 2 K ,ml)(mt) 
+H[' (m 2 K ,m 2 K ,m 2 vl ml)(2ml) 

i tjF' ( 2 2 2 2\fo 4\ 

+H 21 (m w , m„, m^ mj (3m J 
+<(m 2 ,m 2 ,,<;m 2 )(-3/8mt) 

i tjF' I 2 2 2 2\/o 4\ 

+H 21 (m K , m w , m K ; mj (3m J 
+tf£(m 2 ,m 2 m 2 ,;m 2 J(9/8mt) 
+H F (ml, ml, ml; m 2 )(5/3 + 1/18 g 2 ) 
+tf F (m* , m|, m 2 ,; m 2 ) (1/12 m 2 + 1/12 g 2 ) 

+H F {m 2 K , ml, m\; m\) (15/32 m 2 + 5/96 (2m 2 K - m 2 n ) - 5/48 q 2 ) 
+H F (m 2 K , m\, ml; ml) (-5/48 m 2 K ) 
+H F (m 2 K , m\, m^; m 2 ) (-1/2 m 2 - 1/16 q 2 ) 
+H[(ml,ml,ml;ml)(-3ml- 1/3 q 2 ) 
+H[(ml,m 2 K ,m 2 K ;ml)(l/24q 2 ) 

i tjF / 2 2 2 2\f 2 i/o2\ 

+H 1 [m K , m^,m K ; mj (-m w - 1/8 g ) 
+H F {m 2 K , m\, m 2 ; m 2 )(m 2 + 1/8 g 2 ) 
+ H 2 F i( m li m li m l'i m l)(3ml + 1/6 g 2 ) 
+tf£(m 2 , m 2 ,, m 2 ,; m 2 ) (-3/8 m 2 - 1/48 g 2 ) 
+H F 1 (m 2 K , ml, m\; m\) (53/16 m 2 + 1/16 g 2 ) 
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m\, ml; m 2 ) (-5/16 m 2 + 5/48 g 2 ) 

m 3 K , m\; m 2 )(9/8 m 2 n + 1/16 q 2 ) . (A.2) 

Secondly we have the contribution from diagram (b) in Fig. || Most terms in the 
diagram can not be reduced to the previous set of integrals and an independent basis is 
needed. In terms of these we obtain the following result 

cm / 2 2 2 2 22 2\/r / n 4 >r /o 2 2\ 

^yv = Vii{m v ,m v ,m v ,m v ;m v ,q , mJ(5/2 m w - 7/3 g mj 

t; / 2 2 2 2 2 2 2 \ / 4 o /o 2 2 , i /i o 4\ 

+Vu{m n ,m n ,m K ,m K ;rri n ,q ,mj{m n - 2/3 q m w + l/12g ) 
+V n (ml, ml, m 2 , m 2 ; ml, q 2 , m 2 )(l/18 m 4 ) 

+V u {m 2 K , m\, ml, m 2 K ; ml, q 2 , m 2 )(3/2 m 4 n - 17/12 q 2 ml + 1/6 g 4 ) 
+V 11 (m 2 K , m\, m\, m 2 ; ml, q 2 , m 2 )(2/3 m 4 - 2/3 q 2 m 2 n + 1/8 g 4 ) 

T/ / 2 2 2 2 22 2\/ r 2 , 2\ 

+V r 2i(m 2 r ,m 2 r ,m 2 Y ,m 2 r ;m 2 r ,g 2 ,m 2 r )(-2m 2 r + 2/3 g 2 ) 
+V21 (m^ , , ml , m\ ; ml , q 2 , ml)(-4ml + 4/3 q 2 ) 
+V 21 (m 2 K , m 2 K ,m 2 K ,m 2 ;ml,q 2 , m 2 )(-2m 2 + q 2 ) 
+V 2 2(ml,ml,ml,ml;ml,q 2 ,ml)(-6ml + 10/3 g 2 m 2 ) 
+V" 2 2(m 2 ,m 2 r ,m 2 Y ,m 2 r ;m 2 r ,g 2 ,m 2 )(-2m^ + 4/3g 2 m 2 - 1/6 q 4 ) 
+V2 2 (m 2 K ,m 2 K ,ml,m 2 K ;ml,q 2 ,ml)(-4ml + 11/4 g 2 m 2 - 1/3 q 4 ) 
+V 2 2(m 2 K ,m 2 K ,m 2 K ,m 2 ;ml,q 2 ,ml)(-2m 4 + 5/3 g 2 m 2 - 1/4 q 4 ) 

,/ / 2 2 2 2 22 2\/r- /o 2 2 i 1 /o 4\ 

+V 24 {m 7T ,m 7r ,m 7r ,m 7r ;m 7r ,q ,mj(5/3g m w + l/2q) 

ir / 2 2 2 2 2 2 2\/i/o2 2\ 

+V 24 {m n ,m 7r ,m K ,m K ;m 7r ,q ,mj(l/3g mj 
+V" 24 (m 2 r , m 2 K , ml, m 2 K ; m 2 n , q 2 , m 2 )(5/6 g 2 m 2 ) 

,/ / 2 2 2 2 2 2 2\/i/o22\ 

+V 2i {m K ,m K ,m K ,m v ,m n ,q ,mj(l/3g mj 
+V 25 (m 2 n , ml, ml, m 2 ; m 2 , g 2 , m 2 )(-4m 2 + 4/3 g 2 ) 

,/ / 2 2 2 2 22 2\/ o 2 i 1 /o 2\ 

+V 25 {m n ,m 7T ,m K ,m K ;m 7T ,q ,m 7r )(-2m 7r + 1/2 g ) 
+V 25 (m 2 K ,m 2 K ,ml,m 2 K ;ml,q 2 ,ml)(-3ml + 17/12 g 2 ) 

ir / 2 2 2 2 2 2 2 \ / o 2 , 2\ 

+V 25 {m K ,m K ,m K ,m v ,m n ,q ,m^)(-2m^ + q ) 
+V 26 (ml,ml,ml,ml;ml,q 2 ,ml)(-4ml + 7/3 g 2 m 2 + 1/3 g 4 ) 

, T/" / 2 2 2 2 2 2 2\/o4,22\ 

+V 26 {m^,m n ,m K ,m K ;m 7V ,q ,m 7r )(-2m n + q mj 
+V 26 (m 2 K ,m 2 K ,ml,m 2 K ;ml,q 2 ,ml)(-3ml + 13/6 g 2 m 2 - 5/24 g 4 ) 
+V 26 (m 2 K ,m 2 K ,m 2 K ,m 2 -,ml,q 2 ,ml)(-2ml + 3/2 g 2 m 2 - 1/8 g 4 ) 
+V 29 (m 2 n , ml, ml, m 2 ; m 2 , g 2 , m 2 )(4/3 g 4 ) 

if / 2 2 2 2 2 2 2\/i ;< 4\ 

+V 2 g{m 7n m^,m K ,m K ;m 7: ,q ,mj(l/4g ) 
+V 29 (m 2 K , m\, ml, m\; ml, q 2 , m 2 )(7/24 g 4 ) 
+V 29 (m 2 K , m 2 K , m 2 K , m 2 ; m 2 n , q 2 , m 2 )(l/4 g 4 ) 

ir / 2 2 2 2 22 2\/p 2 o 2\ 
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+ ^31 
+ V 31 
+ V 31 
+ V 32 
+ V 32 
+ V 32 
+ ^33 
+ ^33 
+ ^33 
+ ^33 
+ ^35 
+ ^35 
+ ^35 
+ ^35 
+ ^36 
+ ^36 
+ ^36 
+ V 36 

+ ^37 
+ ^37 
+ V 37 

+ ^37 
+ ^38 
+ ^39 
+ ^39 
+ ^39 
+ ^39 

+ ^310 
+ V 3W 
+ V 3W 
+V311 
+ V 311 
+ ^311 
+V311 
+ ^313 



2 2 2 2 2 2 2\/o 2 2\ 

m w , m w , m K , m K ; m w ,q , m J (3m w - g ) 

2 2 2 2 2 2 2\/r 2 o 2\ 

m K ,m K ,m 7T ,m K ;m 7T ,q ,mj(6m x - 2q ) 
m^, m\, m 2 K , m 2 ; ml, q 2 , ml) (9/2 m 2 - 3/2 g 2 ) 

2 2 2 2 22 2\/i/o2\ 

m 7T ,m 7T ,m K ,m K ;m 7T ,q ,mj(-l/3g ) 

2 2 2 2 22 2\/o/o2\ 

m K ,m K ,m n ,m K ;m n ,q ,mJ(-2/3g ) 

2 2 2 2 22 2\/i/r>2\ 

m^m^m^m^m^q ,mj(-l/2g ) 

2 2 2 2 22 2\/r> 4 2 2\ 

m^m^m^m^m^q , mj (2m,, - g raj 

2 2 2 2 2 2 2\/ 4 r> /O 2 2 , i /i r> 4\ 

m^,m^,m K ,m K ;m^,q ,m ff )(m w - 2/3 g ra^ + 1/12 g ) 
m 2 K ,m 2 K ,ml,m 2 K ;ml,q 2 ,ml)(2ml - 4/3 g 2 ra 2 + 1/6 g 4 ) 

2 2 2 2 2 2 2\/o /o 4 2 2 , -i /n 4\ 

ra^ra^ra^ra^m^g ,m 7T ){3/2m 7T - g ra^ + 1/8 g ) 

2 2 2 2 2 2 2 \ / i /<-> 4\ 

m^m^m^m^m^q ,raj(-l/2g ) 

2 2 2 2 2 2 2\/i/o22\ 

m^,m^,m K ,m K ;m^,q ,raj(-l/3g raj 

2 2 2 2 2 2 2\/o/o2 2\ 

m K ,m K ,m^,m K ;m^,q ,raJ(-2/3g raj 

2 2 2 2 2 2 2\/i/o22\ 

m K ,m K ,m K ,m ri ;m w ,q ,raj(-l/2g mj 



7T ' 7T ' 7T ' TV' TV ' 

,2 . 



g 2 X)(-l/2g 4 



2 2 2 2 2 2 2\/i/io4\ 

m^,m^,m K ,m K ;m^,q ,raj(-l/12g ) 

2 2 2 2 22 2\/ -i / c 4\ 

m K ,m K ,m„,m K ;m„,q ,raj(-l/6g ) 

2 2 2 2 22 2\/i/o4\ 

m K ,m K ,m K ,m v ,m n ,q ,mj(-l/8g ) 

2 2 2 2 22 2\//|2\ 

m n ,m 7T ,m 7T ,m 7T ;m 7T ,q ,raj(4raj 

ml, ml, m 2 K , m 2 K ; ml, g 2 , m 2 )(2m 2 — 1/2 g 2 ) 

2 2 2 2 2 2 2\(a 2 2\ 

m K ,m K ,m^,m K ;m^,q ,m 7r )(4m 7r -g ) 

2 2 2 2 2 2 2\/o2 o//i2\ 

m K ,m K ,m K ,m v ;m 7T ,q ,m 7T ){3m 7T — 3/4q ) 

2 2 2 2 2 2 2\/<-> 2\ 

m n ,m^,m^,m^,m^,q ,raj(2g ) 

2 2 2 2 2 2 2\/o 2 a 2\ 

m^,m^,m^,m^;m^,q ,m^)(8m 7r -4g ) 

2 2 2 2 2 2 2\l a 2 o/o2\ 

m^,m^,m K ,m K ;m^,q ,mj(4m ff - 3/2 g ) 

2 2 2 2 2 2 2\/o 2 o 2\ 

m K ,m K ,m^,m K ;m^,q , mj{8m v - 3g ) 

2 2 2 2 2 2 2\/p 2 n / A 2\ 

m K ,m K ,m K ,m r} ;m n ,q ,mj(6m n - 9/4 g ) 

2 2 2 2 2 2 2\/o/o2\ 

m n ,m n ,m n ,m n ;m n ,q ,raJ(-2/3g ) 

2 2 2 2 22 2\ / r /c 2\ 

m K ,m K ,m w ,m K ;m n ,q ,raJ(-5/6g ) 

2 2 2 2 22 2 \ / 1 /o 2\ 

2 2 2 2 2 2 2\//i 4 ,1/0 2 2 i /o 4\ 

m n ,m n ,m n ,m n ;m n ,q ,m 7r )(4m 7r - 4/3 g - 1/3 g ) 

2 2 2 2 22 2\/<->4 2 2\ 

m 7T ,m 7T ,m K ,m K ;m 7T ,q , m J (2m w - q raj 

m\, m\, ml, m 2 K ; ml, q 2 , m 2 )(4m 4 - 29/12 g 2 m 2 + 5/24 g 4 ) 

2 2 2 2 2 2 2\/o 4 ^//i 2 2 , i /n 4\ 

K ,m K ,m K ,m v ;m n ,q , raj (3m w - 7/4 g + 1/8 g ) 



m 



q 2 ,ml)(2q 2 ml - 4/3 g 4 ) 
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+V 313 (ml, ml, m\, m 2 K ; m 2 n , q 2 , m 2 )(-l/4 g 4 ) 
+V 3X3 {m 2 K , m\, ml, m 2 K ; ml, q 2 , m 2 )(-7/24 g 4 ) 
+V 3 i 3 (m 2 K , m\, m\, m\; ml, q 2 , m 2 )(-l/4 g 4 ) 
+V 315 (ml,ml,ml,ml;ml,q 2 ,ml)(-2/3 q 2 m 2 n - 2/3 g 4 ) 
+V 315 (ml, ml, m 2 K , m\; ml, q 2 , m 2 )(-l/4 g 4 ) 

ir / 2 2 2 2 2 2 2\/ r / c 2 2 1/10 4\ 

+V 315 {m K ,m K ,m 7r ,m K ;m 7T ,q ,mJ(-5/6g m w - 1/12 g ) 

/ 2 2 2 2 22 2\/ i /r> 2 2 1 /o 4\ 

+V 315 {m K ,m K ,m K ,m TI ;m 7: ,q ,mj(-l/2g - 1/8 g ) 

ir / 2 2 2 2 2 2 2 \ / r /o 4\ 

+V 316 {m n ,m n ,m n ,m n ,m n ,q ,mJ(-5/3g ) 

T/ / 2 2 2 2 22 2\/ i /r> 4\ 

+K 3 i 6 (m 7r ,m 7r ,m^,m K ;m^,g ,m 7r )(-l/2g ) 
+V 3W (m 2 K , m\, ml, m 2 K ; ml, q 2 , m 2 )(-7/12 g 4 ) 

ir / 2 2 2 2 22 2\/ i in 4\ 

+V 31 e(m K ,m K ,m K ,m v ;m 7r ,q ,mj(-l/2g ) 

iTA/2 2 2 2 2 2 2\//i2 r>2\ 

+V3i7(w 7r ,m 7r ,m 7r ,TO 7r ;m 7r ,g , m 7r )(4m 7r - 2g ) 

ir / 2 2 2 2 22 2\/o 2 2\ 

+^ 3 i7(m 7r ,m 7r ,m^,m x ;m^,g ,m 7r )(2m 7T -q ) 
+V 3 i 7 (m 2 r ,m 2 Y ,m 2 ,m 2 r ;m 2 ,g 2 ,m 2 J(3/2m 2 - 3/4 g 2 ) 

T/ / 2 2 2 2 2 2 2\/o2 o 2\ 

+V3i9(TO 7r , m^m^m^ ;m„,q ,m n ){8m 7r -2q ) 
+^(3i9(m 2 , m^m^m^m^g^m^m 2 - g 2 ) 

ir / 2 2 2 2 22 2\/o 2 o ia 2\ 

+V 319 {m K ,m K ,m 7r ,m K ;m n ,q , raj^m^ - 3/4 g ) 

T/ / 2 2 2 2 22 2\/o 2 o //i 2\ 

+V 319 {m K ,m K ,m K ,m rj ;m w ,q ,m n )(3m n - 3/4 g ) 

ir / 2 2 2 2 22 2\/ / i 4 <-, 2 2 \ 

+^32i(m 7r ,m 7r ,m 7r ,m 7r ;m 7r ,g ,m n )(Am 7r -2q mj 

ir / 2 2 2 2 2 2 2\/r>4 2 2 \ 

+V 32 i{m w ,m n ,m K ,m K ;m 7T ,q ,m 7r )(2m 7r -g mj 
+V32i(m 2 Y,^x,m 2 ,^x;^ 2 ,g 2 ,^ 2 )(3/2m 4 - 3/4 g 2 m 2 ) 
+V 32 i(m 2 K ,m 2 K ,m 2 K ,m 2 i ;ml,q 2 ,ml)(3/2ml - 3/4 g 2 m 2 ) 

iTA/2 2 2 2 2 2 2 \ / 4\ 

+V 323 {m 7T ,m 7r ,m 7r ,m 7r ;m 7r ,q ,mj(-g ) 

ir / 2 2 2 2 2 2 2w 1/0 4\ 

+V 323 {m n ,m n ,m K ,m K ;m 7r ,q ,mj(-l/2g ) 
+^32sOL m 2 K , ml, m 2 K ; m 2 , q 2 , m 2 )(-3/8 g 4 ) 
+ V323 (m 2 K , m 2 K ,m 2 K ,m\;ml,q 2 , m 2 ) (-3/ 8 g 4 ) 

IT//2 2 2 2 2 2 2 \ / 4\ 

+V 325 {m 7T ,m 7T ,m 7T ,m 7T ;m n ,q ,mj(-q ) 
+V 3 2 5 (ml , ml, m\, m 2 K ; ml, q 2 , m 2 )(-l/2 g 4 ) 

ir / 2 2 2 2 2 2 2\/o/ 4\ 

+V 325 {m K ,m K ,m 7r ,m K ;m n ,q ,mJ(-3/8g ) 

+V 3 25(m 2 K ,m 2 K ,m 2 K ,m 2 ;ml,q 2 ,ml)(-3/8q 4 ) , (A.3) 
where the notation can be read from App. 0. 
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B One-loop integrals 

Through the calculation one has to use one-loop integrals of one, two and three point 
functions. The latter disappear after mass renormalization and the use of some recursion 
relations. We only have to deal with the following set of functions - in the remainder we 
use d = 4 — 2e. 



A(m\) 
B(ml,m 2 2 ,p 2 ) 
B^(m\,m\,p 2 ) 

B^(m\,m\,p 2 ) 

2 „2\ 



1 f d d q 



i J (27r) d q 2 — mf ' 
1 f d d q 1 



i J (27r) d (q 2 — m\)((q — p) 2 — ml) ' 
d d q q^ 



i J (27r) d (q 2 — m\)((q — p) 2 — ml) 
p^B^m^m^p 2 ) , 

d d q q^q v 



B lxva (m 1 ,m 2 ,p 



i J (27i) d (q 2 — ml)((q — p) 2 — ml) 
Vv>VvB 2 \ (ml ,m 2 2 ,p 2 ) + g flv B 22 (m\,m\,p 2 ) 
d d q q^q v q a 



ij (2n) d (q 2 — m\)((q — p) 2 — m 2 ) 
PvPvP a B 31 (ml,ml,p 2 ) + (p^g va +p v g m + p a g^)B 3 2(ml,ml,p 2 ) . 

(b.i; 



An expansion in e leads to the following series 



,2 



m 

16tt 2 

.2 2 



A (™1) = T7 ^Xo + A(m 2 1 ) + eA t (m 2 1 ) + 

ITT 
1 



Bi^m^m^Sn) = Y^pole^ + B ij (m 1 ,m 2 , s w ) + eB ij (m 1 ,m 2 ,s^) + . . . , (B.2) 

with A,Bij defining finite quantities and where "pole^" denotes the singular part of each 
of the Bij functions, 

pole = A , pole 1 = y, pole 21 = y , pole 22 = y (m\ + m\ - y ) , 

P° le 3i = Vo\e 32 = ^(2m 1 + Ami - s w ) , ( B - 3 ) 

with 

A = - + ln(47r) + 1 -7. (B.4) 
e 

After some simpler algebraic manipulation, the functions defined in Eq. ( |B.1| ) can be 
related to the basic integrals A(m 2 ) and Bi(m 2 ,ml,p 2 ) through the identities 
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B 31 (ml, ml s n ) = ±- (AK) - (m 2 -m\- s^)B 2l {ml < **) - ^K, ml , 

I / 777^ 777^ \ 

B^m^mls^) = —^-—A(ml) + ^A(ml)-(ml-ml-s^)B 2 2(ml,mls^)J, 
B 21 (ml,mls«) = -(A(ml)+mlB(ml,ml,s^-dB 22 (ml,mls n )), 

-(sv + ml-mDB^mlmlsj), 
B x (m\,mls^) = — (A(m 2 2 ) - A(m\) + {m\ — m\ + s ff )S(m?, m^O) > 
fl(m?,m?,0) = ^A(m]). (B.5) 

Notice that the previous identities are only used at the final numerical level in order 
to avoid cancellations between different terms occurring in the form factors. The explicit 
expressions are 

A{m{) = -j^mK), 



B[m x ,m 2 ,v 



1 m\ ln(mf) — m\ ln(m|) 
167r 2 m\ — m\ 

1 / /A £\ , m? i/ , Qo 2 + z/) 2 - A 2 . 



(32vr 2 ) V V P 2 A,/ m 2 p 2 (p 2 - z/) 2 - A 2 

with A = m\ — m 2 , £ = m\ + m 2 , and z/ 2 = [p 2 — (mi + m2) 2 ][p 2 — (mi — TO2) 2 ]. Similarly 
combining Eqs. (|B.2| ) and ( p.5|) one can obtain the expressions of the e terms in the 
expansion as functions of A and B . A straight forward calculation leads to 

l^A\m\) = m 2 [^ + i + ^ + iln 2 (m 2 )-Cln(m 2 )], 

_ q2 i ^2 _ i rl 

16n 2 B\ml,ml,s n ) = —-- + — + (C - l)5(m 2 , m 2 , s,) + - / dx ln 2 (m 2 ) , 

(B.7) 

with C = ln(47r) + 1 — 7 and m 2 = xm\ + (1 — x)m\ — x(l — x)s 7T . The definition ( |B.1| ) 
is /j, independent. The subtraction procedure will always provide the scale /1 in the correct 
fashion to give all logarithms dimensionless arguments of the type ln(m 2 //i 2 ). 
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C Vertex integrals 



In this appendix we display the vertex type integrals. They were needed to calculate the 
form factors |], but we never showed their definition nor how they appear in the 
expressions. The pion vector form factor has a manageable expression, that allow us to 
show them. It thus makes sense to display their definition and evaluation here. 
We define 

1 f d d r d d s X 



i 2 J (2n) d (2n) d (r 2 -m 2 )((r-q) 2 -m 2 )(s 2 -m 2 )((r + s-p) 2 -m 2 ) ' ^ 

In the calculation one encounters integrals of this type with up to three integrated 
momenta in the numerator. We Lorentz decompose them into scalar functions. All these 
functions, referred to as vertex integrals, depend on m 2 , m|, 777J, m 2 , p 2 , q 2 , (p — q) 2 , ji 2 . The 
scale [L appears always in connection with the //-independent integrals as defined in ( |U.1| ) 
to make all arguments of logarithms dimensionless. In the following Lorentz decomposition 
the dependence on these arguments inside the Vi functions has been suppressed for brevity. 



yr 

I" 

V s 
>' 

yrr 
yrs 

[IV 



V = ((1)), 

((Vk)) = 9fiuV2i + PnPvV 2 2 + q^q v V23 + {p^qu + q^Pu)V 2 i , 
: ((v^)) = g^V 25 + p^p v V 2& + q^qyV 27 + q^p v V 2 % + p^q v V 29 , 
v £ = ((V*)) = 9ftuVa.o + P&uVm + q^uV 212 + (q^p u + p^q v )V 213 , 
v iZa = ((v^)) = (g^p a + g m p v + g va p^)V 3 x + (g^q a + g m q v + fiwg M )V 32 

+PpPuPaV 33 + q^qaVu + (PpPuQa + PpQvPa + q^PuPa)V 3 5 

+(QnQuPa + q^q a + p^q a )v 36 , 
v JZa = ((r^r u s a )) = g^PaV 37 + g^q a V 3& + {g m p v + g va Pn)V39 

+q l j,quPaV3u + {p^q u + q^Pu)p a V 3 i 5 + {p^q v + q^p v )q a V 31& , 

= (( r f* S " s a)) = Pn9vaVai7 + ? M flW^318 + (g^uPa + g^aPu)V 319 

+ (g^q a + g m qu)V 32 o + P^P a V 32 i + q^quq a V 322 + Pf,q„q a V 323 
+q^PuPaV 32 4 + P^{Puq a + quPa)V 325 + q^{p u q a + quPa)V 32& 



(C.2) 



This set of 44 functions is not obviously a minimum set and some relations can be 
found between then. Instead of reducing to a basic set, just by simply manipulations like 
contracting with the external momenta or a g^ tensor, we use those as cross-check of our 
numerical integrals. In addition one can relate functions with different arguments with 
some redefinitions of momenta. For instance the substitution r — > —r' + q; s —>■ —s' relates 
functions with the canonical arguments m\, rnj, m\, m 2 ,p 2 , q 2 , (p — q) 2 ,^ 2 to those with 
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arguments m|, m\ , m'g, m 



4 , (q — p) 2 , q 2 ,p 2 , /i 2 , while r — > r'; s — > —r' — s' + p relates them to 
those with arguments m 2 , mj, mf, m|, p 2 , g 2 , (p — g) 2 , /i 2 . As an especial case we find a huge 



simplification in our results when we deal with the substitution r 



r ; s 



-r — s + p 



for the more restrictive case m 2 
various integrals in Eq. ( JCJ.2| ) 



obtaining the following set of relations between the 



V 



yrs 



;(PuV:-V 



rr \ 



jrrrs 



( n v rr - V 



rrr ^ 



(C.3) 



In order to evaluate these functions we used the methods developed in 27|, where 
we refer the reader for a more extensive treatment. Here we only repeat the basic steps 
- our functions differ slightly from the ones in |27| since we stay in Minkowski space 
throughout. 

We first combine the first two propagators with a Feynman integration in Eq. ( |C . 1| ) 



«*» = /V 1 ! iir dia 



z 2 J (2tcY (2nY [(r - yq) 2 - ((1 - y)m\ + ym\ - y{\ - y)g 2 )] 2 

' (C4) 



(s 2 — m|)((r — yq + s — (p — yq)) 2 — m 2 ) ' 
and we shift the r integration to r — yq. All remaining integrals can be written in terms of 
2 2 2 ,2a 1 [ d d r d d s (r ■ fc)*(s • jty 



,m 3 ,m 4 ,A;)- . 2 J [ r 2 - m 2 ) 2 {s 2 - ml){{r + s + k) 2 - m 2 ) (C ' 5) 

where in our case ij is reduced to the values 00,10,01,20,11,02,30,21 and 12. We also 
have defined 

m 2 = (1 — y)m\ + ym\ — y(l — y)q 2 , k = -p + yq . (C.6) 

The infinite parts of these functions can be obtained analytically and their finite parts can 
be written in terms of 9 functions, hi, that have a fairly simple integral representation - 
see |?7] for an explicit representation of these functions. So in the end all the functions 
appearing in the vertex diagrams are obtained as a double integral. The singularities in 
these integrals can be avoided by deforming the integration paths , this also serves as a 
check on the numerical code. 

Following the previous steps is straight forward to obtain the relations between the 
vertex, Vi, and the basic integrals, P^ c , in terms of the variables defined in Eq. ( |U.6|) - in 
what follows and for simplicity Vi should be read Vi(m 2 , m 2 , 7??J, m 2 ;p 2 , q 2 , (p — q) 2 )- 



V = {P™} [m 2 ,ml,ml;k 



piu 

Vn = \~^r \ [m 2 ,m 2 3 ,ml;k 2 ] 



2. 1.21 

y 

10 



y 

plO 

-To 



VL. = {y(PZ-^)}y,ml,ml;k 2 \ 
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V* = {-^[m 2 ,m 2 ,m 2 ;A; 2 ], 
r P 01 -i 

V " = {y^\jm 2 ,mlml;k 2 }, 

f i p20 p20 . 

, „2 p20 p20 plO ,. 

V 25 = ^^{-2P^-P^-2P^-(e + m ^ + m l- m l)PZ 
+B(m 2 , m\ 0)(A(mj) - A(m 2 )) - ^)) [m\ m§, m 2 ; A; 2 ] , 

Vfc = { ^ _ 1} + fin + 2P°\ + (A: 2 + m 2 + m 2 - m 2 )^ 

-P(m 2 , m 2 , 0)(A(m 2 ) - A(m 2 )) + 2^) + ^}jm 2 , m 2 , m 2 ; A; 2 ] , 

^ = { 2fc2( f _ 1} (2P& + Pm + 2P°\ + (A; 2 + m 2 + m 2 - m 2 )^ 

-P(m 2 , m 2 , 0)(A(m 2 ) - A(m 2 )) + 2^) + ^(P^ + ^)}^ 2 , m 2 , m 2 ; A; 2 ] , 

V» = { 2k2( l _ 1} ( - 2P, 1 " - C - 2P°\ - (A; 2 + m 2 + m 2 - m 2 )^ 

+P(m 2 ,m 2 ,0)(A(m 2 ) - A(m 2 )) - 2^) - |(P 2 A + ^)} Jm 2 , m 2 , m 2 ; A; 2 ] , 
V 29 = { 2fc2( ^ _ i} ( - 2P 2 \° - PZ - 2P°\ - (A; 2 + m 2 + m 2 - m°)i*° 

+P(m 2 ,m 2 ,0)(A(m 2 ) - A(m 2 )) - 2^) - ^P 2 \\}^[m 2 , m 2 , m 2 ; A; 2 ] , 
V 2W = ^^{m 2 P 2 7 1 + A(m 2 )P^^ 

^2ii = {-^^(-m 2 P 2 °{\^^ 

^212 = {^^(-m^i^^^ 
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Vs2 = (^{y(™ 2 ^i + ^ 

, p30 q p30 . 

v » = { - + fi^tt) < m2 ^» + p ''° - ¥ ] )y< < < e] ' 

( / 3 3 p 30 3 

T/ L.3/ pOO I 71)10 I p20 , £_211 , /- ^pOO _ pOO 

^34 - (y (P 211 + ¥ P 2 U + F ^211 + + fc 2 (n " Pill 

+^(-^ 2 ~ Pin + 4°i + ATi))) ml ml k 2 } , 

+ l(-3m 2 P 2 \° 1 - 3< + < + 3^))) }^[m 2 , m 2 , < k 2 } , 
( v 2 / 2 P 30 1 

T/ _ J « / plO z p20 ^211 , 1 /q 2 pOO , 9 pOO 

+^(3m 2 P 2 \° 1 + 3P& - 2P 2 2 ° - 3^))) } y [m 2 , m 2 , m 2 ; A; 2 ] , 

Vss = {^^yP^ + P^-^y^lmlk 2 }, 

V39 = { 2k\n - 1) ( (fc2 + m2 + m l~ p 2ii + Plii + 2P 2 \\ + 2P 2 2 ° 
+2^)}jm 2 ,m 2 ,m 2 ;A; 2 ], 

-PS - 2P 2 °A + l(-(m 2 + ml + 3k 2 - m 2 )P 2 \° - P™ - 4P 2 \\ - 2P 2 2 1 ° 1 
-2^))} y [m 2 ,mlmlk 2 ], 
Vsn = {^(-^ + j^yP^ + Pn\-(k 2 + m 2 + ml-m 2 4 )P^ 
-2P 2 \\ - 2P 2 2 1 ° 1 - 3^)) } y [m 2 , ml ml k 2 } , 
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= {^(<+2^ + ^ + ^(5(- 2 ,- 2 ,o)(ak)-a ( ^)) 

+ (k 2 + m 2 + m 2 3 - ml)P™ + P™ + 2P 2 °A - ^P°\ - ^ 
+ i_(( 3A; 2 + m 2 + m 2_ ro 2 )p io + pio + 4 pn + 2P 20 

+3^)))}jm 2 ,m 2 ,m 2 ^ 2 ], 
^313 = {|f(^+(^(-m 2 P 2 i-^ 

+2P 2 \\ + 2J*° + 3^|)) } [m\ ml m 2 ; k 2 } , 

= { y ^(-p^-2^-^ + ^(B( m i m 2 MA( m t)-A( m i)) 

-(A; 2 + m 2 + ml - m 2 )^ - P™ - 2P 2 °A + -^(m 2 ^ + P{\\ 

-(3A; 2 + m 2 + m\ - m 2 )^ - - 4P 2 \\ - 2P% X - 3^))) } \m\ m 2 , m 2 ; k 2 } , 

+ \(k 2 + m 2 + mt- ml)P™ + \p^ + P°\ + ^(-m 2 ^ - P°\ 

+ (2A; 2 + m 2 + m\ - m 2 4 )P^ + + 3P 2 \\ + 2P 20 ! + 3^))) }jm 2 , m 2 , m 2 ; A; 2 ] , 

-i(A; 2 + m 2 + m 2 - m 2 )P 2 °{\ - ± J*° - P°\ + ^(m 2 P 2 °A + P°\ 

-(2A; 2 + m 2 + m\ - m 2 )^ - < - 3P 2 \\ - 2P™ - 3^))) }^[m 2 , m 2 , m 2 ; A; 2 ] , 

^sir = j^{l 2 (- m lP^ + ^)}yM,<k 2 }, 

V 318 = -^^{y (A(ml)B(m 2 , m 2 , 0) + m 2 ^ + ^(-^Ti + mfon 

-^ ) )}, [m2 ' m '' m ^ 2] ' 
V 319 = l _ - [ - A(m 2 4 )B(m 2 , m 2 , 0) + 1 ((A: 2 + m 2 + m 2 - m 2 ^ 

+P 1 ° 1 \ + 2i** + 2P 2 \\ + 2^) } Jm 2 , m 2 , m 2 ; A; 2 ] , 
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^20 = 2(w 1 _ 1} {y(A{ml)B{m\ m\ 0) + ^(-(k 2 + m 2 + mj - m 2 )P 2 n 

-P n \ - 2P 2 ° n - 2P 2 \\ - 2^)) }^[m 2 , mi m 2 ; k 2 } , 
^32i = {^(-^ + ^^(A(mDP(m^m^0) + l(-(A; 2 + m 2 + m 2 -mDP 

-An - 2P 2 Ti + - 2P 2 \\ - 3^))) } Jm 2 , m 2 , m 2 ; A; 2 ] , 

^322 = {|J(^ + ^ + ^(-2^P(^^0)-m 2 P- 

+ 1((A; 2 + m 2 + m\ - m 2 )P 2 ° n + P lTl + 3P 2 ° n - m 2 ^ + 2P 2 \\ 



where 



211 



+3^)))}jm 2 ,m 2 ,m 2 ;A; 2 ], 
V-323 = {|J(-^+(^(^ 2 )P(m 2 ,m 2 ,0) + l(-(A; 2 + m 2 + m 2 -m 
-An - 2P 2 Ti + " 2P 2 \\ - 3^))) } [m 2 , mj, mj; k 2 } , 

+ i_( (A; 2 + m 2 + m 2 _ ^pcn + pea + g^H _ m 2p 2 io + 2p n 
+3^)))}jm 2 ,m 2 ,m 2 ;A; 2 ], 



pl2 



-P n \ + 2P 2 ° n - mjP™ + 2P 2 \\ + 3^))) } [m 2 , m 2 , m 2 ; k 2 } , 



A; 2 

V» = {g(-^-^ + ^(2A(^(^^0) + m 2 P 2 n 

+ 1(-(A; 2 + m 2 + mj - m 2 )P 2 ° n - P n \ - 3P 2 ° n + mjP^ - 2P 2 \\ 



" 3 ^ )) )}, [m2 ' m '' m ' ;A;2] ' (C - 7) 



{•••!, • I' dy{...}, 
Jo 



(C.8) 



and we have display the P l { arguments outside the brackets. 



abc 

n b( 

reason of choosing Pgii as a fraszs instead of the simpler P^ is the infra-red behaviour in 



The functions P^ can be obtained in terms of the P^u via recursion relations. The 



43 



the last set. As is evident the Pi U have a direct relation with the Hi functions defined in 
@, and we have checked that they agree with each other. The relations between the P^ m 
functions with the 9 hi is given in Minkowsky space-time by (see [27]) 



A 5 



(16^)^(1, 2, 3; fc 2 ) = 'f(ml + mj + ml) 



+Ai 



k 2 1 1 1 

"J + mli- - hii) + m 2 (- - ln 2 ) + m 2 (- - ln 3 ; 



/ 9 1 \ 1 7T 2 

+ k 2 ^ - - + - ln X -/ii ; i23 + /i2;123j + m l ( 2 + + ln l ~ 

1 7T 2 1 7T 2 

+m % + 12 + 1112 ~ 1112 +hl ' 213 " > + m ^2 + 12 + ^ ~ 1113 +/ll;312 ' 1 



A 



(16^)^(1, 2, 3; A; 2 ) = -'^k\m{ + m 



_^Ifc 2 ( - — + m?(- - lm) + ^1 + m 2 (- - ln 3 ) 
2 V 12 1V 4 ' 2 4 y 

-fc 4 ( - g + g ln 2 --/il;213 + g^2;213 - j^^is) 
,2/ 2/ 3 ^ 1, 2 1, 1, 2, 

-fc mf ( 1 h - lnf — lni +-^1-123 + ~h 

V u 16 24 2 1 4 6 ' 3 

7 1 2 
_m 2(g - 2 1112 + 3 ^ 1;213 ~~ /l2 = 213 )) 

2/ 3 71-2 1, 2 1, 1, 2, 
- m 3(-^ + 24 + 2 3 ~4 3 + 6 1:321 + 3 3;321 ' 



^(1,2, 3; k 2 ) =^(2,1,3;^), 
(levr 2 ) 2 ^^!^^;^ 

(16^)^(1,2,3;^ 



^A 2 - ^A x (i + lm) - i + ^ + in 2 + im +/ii ; i 23 , 

--^k 2 + k 2 ( - j| + ^ lni -/ii ; i23 + ^123) , 



^2 + p Ai( 3 + l hli) 



A; 2 /9 7T 



(leTT 2 ) 2 /*^!,^;* 2 



12 

^2,2/ 2,1 2,1 2\ 

— fc (m x + -m 2 + -m 3 ) 



+Ai^ ( ^(7 " lni) + \m 2 2 ( 3 - - ln 2 ) + \m\{\ - ln 3 )) 



2/ <mi 1 
2 l 4 



,4/11, 5, 3, 

+ k [ - q + 2 hl , 123 ~ 4 /l 2;123 + 7^4;1: 

)2 /m 2 .l 7r 2 2 1 1 , m 2 .7 7r 2 ,2 s ,, 
+* (T<8 + II + " 2 '"'> + T<8 + 12 + "2 lnj) 

m 3 .7 7T 2 , o 3 . . In, In, In, \ 

+ ^(g + ^ + ln 3 - 3 ll13 ) + ^l^ 123 + Y^^l^lS + 4^3^312 J 
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(16^)^(1, 2, 3; e) = -^k\m\ + mf) 



+ ^(^ + m ?(l ni -i)-!=I +m I(.n3-^)) 
.x/19 1, 1, 3, , 3, 



1 3 '6 \ 

+^(43 - ^ lni +-/ll;123 - g^2;123 + ^3;123 - ^^5;123j 
.2 _2 



, 2 / m\ .1 7T 2 , 2 1 , , 

+ I " ~4^8 + 12 + 1 ~2 1 + 1;123) 
m 2 111 1 1 

+ X (4 - 4 ln ? +2 lR l ln 2 -4 ln 2 + ln 2 - g (^l;123 + ^l;21 3 )) 
m 3 / 9 71-2 1 1 2 1 , 3 1 2 1 , , x\ 

+~§~( _ 4 ~ y + 2 1 ~ 1 3 ~2 3 + +/l1 ' 123 ~ /ll ; 312 ' ) J > 

(16^)^(1, 2, 3; A; 2 ) = ^ + ^ + m 2 ) - ^A; 2 (^(ii + In,) + m 2 /- + In,) 

2, 1 , , ^ A; 4 / 37 7T 2 , 2 11 oL \ 
+m 2 (- + lri! + ln 3 ) j + - ( - — + — + In 2 + — lm +3/* 6;123 J 



k 2 



1 1 + ^ + ln 2 +iln 1 +/i 1;123 ^ 



+m ' ( I + T + ^ ln ' + lni 1113 + ^ lni + \ 1113 + ^ 1113 



(16 7 r 2 ) 2 P 2 T 1 (l ! 2,3;A; 2 ) = -^A;V^ + ^) 



"g^l;123 + "^T^2;123 — 2^ 4;123 ^ 7 ; 123 ) 



(l2 " m * + m2(hl2 "T^ + m ' (ln3 "ilO 

A; 4 / 2/ 29 1 , 1 11 5 

+ yr ( 96-i2 lni -^ i23 + ^ 2 
5 5 

+m 2 (-- + lni -g(/ii ; i23 - h 2 -i2z)) 

2/ 17 ^ 2 1 , 2 5 , 1., 2, 
+m2( 144 " 24 " 2 1112 + 12 1112 "6 (/ll;213 + 3* 3;213)) 

g^3;312)) , 



2/ 17 7T 2 1 ]2 5, 1, 

+m3( 144 "24 "2 1113 + 12 1113 -6^ 312 



(16^)^(1, 2, 3; A; 2 ) = -^m\ + ±k\ - ^ + m\{\ + In,) - ^) 

5 , 1 , 1, 1, 3 , 7 , 

+ H 128 + 48 1111 "12 1112 ~4^ 123 + 3^ 213 + 8 ^ 123 " 12^ 213 
l i ! L 3 , 1, \ ,4/ 2/37 7T 2 1, 2 

-2^123 + 4^4;213 + ^5;123 ~ ^jlM J + k ^ ( — - — - - 1^ 

_ 8~ lni + \^ hl ' 123 ~ 5/l 3;i23)) + ml(- ^ + In 2 ln 2 +^ In 2 

-24^i;2i3 + 3^2;2i3 + g^i;i2 3 ) + m 3 {— - — l ni +- lri! ln 3 -—In 
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+ g m 3 — g"l;123 + 2^^1;312 — 2%;321, 

(167r 2 ) 2 P 2 \ 2 1 (l,2,3;A; 2 ) = ^A; 4 K + m^) 

+ T A;4 ( ~ T + m ' ( 12 " lni) + m ' ( Ti ~ 1113 

, 6 / 43 7 1 1 1 1 7 

+ H " 192 + 72 Ml + 36 1112 ~ 12 15123 " 9 1;213 + 8 2;123 + 36 2;213 

k 1 , !, 3 , 1. \ ,4/ 2/ 31 ^ 

~6 3;123 ~ 24 4;213 + 4 5;123 ~ 4 6:123 + 2 ;123 / V m ^~288 + 48 

1 , 2 1 , 1, 5 , . m 2 1 1 , 2 1 , , 

+ - n, liii +— /ii-123 H ^3-123 ) H 1 — hii — hii mo 

4 1 24 9 ' 18 6 v 8 8 1 4 

+ lni +- In 2 - ln 2 --h 1]123 + —h 1;2 i3 + h 2 - 123 - ^^213 + ^3 ; i23) 
o 4 12 6 

2, 77 ^ 1 , 2 1 , , 3 . 2 1 . 

+mi( 1 lnf +- lni ln 3 H In, ln 3 

3V 288 48 16 1 8 16 3 24 

— — ^l;123 + ^T^l;312 + 7T^3;321 ~ 7^3;123) ) • (C-9) 
24 72 9 6 / 

Where for sake of brevity the following notation should be understood 

PUa,b,c,e) = CK^m^ 2 ), 
hi-abc = hi(ml,ml,m 2 c ,k 2 ) , 

9 

Tfl 

ln„ = ln(-f), 
A 2 = A 2 + (ln(4vr) + 1 - 7 ) 2 , 

Ai = A + 1h(4tt) + 1 - 7 , (CIO) 

and Ao is given in Eq. (p.4j ). We have also employed some direct relations between the 
integral representation of the hi functions 

hl ]a bc — hi. ac b , and h 3 - 123 + h 3 . 132 = ^2;123 • (C-ll) 
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